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THE  EFFECT  OF  AERODYNAILCC  NON-LINEARITIES  AND  CROSS -COUPLINGS  ON  THE 
PERFORMANCE  OF  A  CRUCIFORM  MISSILE 
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W.S.  Brown,  V^.E.  Dean,  P.  Hampton, 
H.  Lewis  and  D.I.  Paddison 


SUIiMARY 


This  Note  describes  an  investigation,  made  with  the  aid  of  TRIDAC  and 
MERCURY  computers,  into  the  effect  of  aerodynamic  non-linearities  and  cross¬ 
couplings  on  the  performance  of  a  cruciform  guided  missile,  including  tho 
effects  of  such  aerodynamic  characteristics  on  stability  and  homing.  Aero¬ 
dynamic  data,  obtained  from  wind  tiannol  tests  of  tho  English  Electric  D.A 
configuration  used  to  develop  TPIUNDERBIRD,  and  of  W,2  (V.R.725)  together  vath 
control  and  guidance  system  data,  wore  employed  to  construct  tho  mathematical 
model . 


The  main  conclusions  reached  were  that,  in  spite  of  some  of  tho  awkward 
aerodynamic  characteristics  assumed,  a  cruciform  missile  provided  with  a 
suitable  feedback  control  system  can  function  satisfactorily  up  to  largo  body 
incidences  (30°)  and  need  not  become  unstable.  The  main  difficulty  posod  by 
the  aerodynamic  characteristics  is  with  respoot to  rolling  behaviour,  in  that 
they  cause  the  missile  to  have  preferred  planes  of  incidence.  IVith  roll  rate 
stabilisation,  this  may  result  in  considerable  rolling  during  homing,  with 
some  increase  in  miss  distance.  A  slow  rate  of  roll,  demanded  continuously 
to  smooth  out  the  effects  of  fortuitous  bias,  also  tended  to  increase  miss 
distance  when  applied  to  a  missile  without  bias.  There  was  some  indication 
that  roll  position  stabilisation  leads  to  mainly  smaller  miss  distances  than 
does  roll  rate  stabilisation  in  homing  on  targets  turning  off  course  at  a 
steady  rate.  The  importance  of  accurately  representing  the  aerodynamics  of 
the  missile  is  emphasised. 

Although  TRIDAC  proved  to  be  capable  of  generating  the  response  and 
stability  characteristics  of  the  missile  accurately,  it  was  not  able  to 
generate  miss  distances  with  tlie  precision  required  to  assess  non-linear 
effects  accurately,  although  trends  and  orders  of  magnitude  wore  correctly 
indicated.  Shortcomings  of  tho  equipment  were  responsible  for  small  enrors 
present  in  tho  axis  transformations  to  which  the  values  of  miss  distance  were 
sensitive.  Consequently,  most  of  tho  results  on  homing  were  obtained 
digitally  on  IffiRCURY. 
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Technical  Note  No#  VifE*5 


1  INTRODUCTION 

The  work  described  in  this  Note  grev;  out  of  a  suggestion  that 
Cartesian  missiles,  required  to  operate  at  high  altitudes,  might  be  in 
trouble  because  of  non-linearities  and  cross-couplings  in  their  aerodynamics* 

^  This  impression  had  gained  currency  following  a  report*^  that  some  American 

missiles  of  cruciform  design  had  performed  erratically  during  their  early 
flight  trials,  having  shown  a  tendency  to  roll  instability  at  large 
incidences,  a  feature  which  was  later  ascribed  to  aerodynamic  characteristics 
/  of  the  type  mentioned  above.  It  was  suggested  that  tho  non-linearities  and 

cross-couplings,  which  had  occurred  in  fixed  wing  missiles  with  rearward 
control  surfaces,  affected  the  pitching  and  yavring  moments  and  could  result 
in  loss  of  stability  in  pitch  or  yaw,  while  the  cross-coupling  also  produced 
a  rolling  moment  which,  under  conditions  of  combined  pitch  and  yaw,  increased 
with  incidence.  The  worst  effects,  however,  were  stated  to  be  confined  to 
incidences  exceeding  15°  and  were,  therefore,  unlikely  to  affect  performance 
at  low  altitudes. 

The  non-linearities  were  attributed  in  the  main  to  partial  shielding 
of  the  control  surfaces  by  the  body,  while  the  cross-coupling  was  thought  to 
be  due,  in  part,  to  vortices  shed  from  the  nose  and  forward  part  of  the  body 
of  the  missile.  Whatever  the  cause,  it  seemed  more  than  likely  that  British 
missiles  of  this  type  would  show  similar  characteristics  and  that  THUrroERBIRD, 
in  particular,  would  be  affected,  since  it  had  been  designed  with  high- 
altitudc  capabilities  in  mind. 

It  v/as  agreed,  therefore,  that  there  should  be  a  general  investigation 
of  the  effect  of  aerodynamic  non-linearities  and  cross-couplings  on  cruciform 
missile  performance,  and  that  a  mathematical  model  of  a  cruciform  missile 
should  be  constructed  on  TRIDAC.  Depending  on  the  outcome  of  the  work,  a 
programiTie  of  flight  tests  was  also  envisaged. 

A  prime  requirement  for  the  simulator  work  was  obviously  the  provision 
of  detailed  information  on  the  aerodynamic  characteristics  of  some  missile. 

In  the  event,  the  only  suitable  data  related  to  the  D.4  test  vehicle  used  to 
develop  THUNDERBIRD,  a  l/l5th  scale  model  of  which  had  been  tested  at  R.A.E. 
This  v;as  not  inappropriate  however,  for  reasons  given  above. 

2 

The  wind  tunnel  tests  had  boon  made  at  Mach  numbers  of  1*58  and  2*02, 
up  to  a  total  incidence  of  22°  at  the  lower  Mach  number,  and  25°  at  the 
higher.  The  control  surfaces  had  been  deflected  singly  and  in  combination, 
during  the  tests,  and  the  model  had  been  rolled  by  stages  through  360°.  A 
fair  amount  of  information  on  cross-couplings  was,  therefore,  available. 

^  This  information  formed  the  basis  of  a  three-dimensional  mathematical 

model  of  the  D.l^.  test  vehicle  which  was  constructed  on  TRIDAC.  To  this  model 
there  was  added,  first,  the  control  system  of  W.1  (THUNDERBIRD)  and,  later, 
that  of  Y/.2  (V.R.725)-  These  hybrid  models  were  used  to  study  the  effects  of 
the  aerodynamic  characteristics  on  stability  and,  in  parallel  with  this  work, 
a  mathematical  analysis^  was  made  of  the  static  stability  characteristics  by 
perturbation  methods.  Later,  it  was  deemed  essential  to  obtain  more 
extensive  aerodynamic  data.  The  development  of  THUNDERBIRD  into  W.2  (V.R.725) 
v/as  in  hand,  and  as  this  involved  certain  changes  in  the  aerodynamic  con¬ 
figuration  of  the  missile,  it  was  clearly  desirable  to  have  up  to  date  data 
for  the  mathematical  model.  A  series  of  wind  tunnel  tests  was  therefore 
arranged  in  the  A.R.A.  supersonic  tunnel  at  Bedford.  These  covered  a  much 
wider  range  of  incidence  and  Mach  number  than  the  D.4  tests  and  included  a 
more  detailed  investigation  of  the  control  surface  characteristics.  The  data 
so  obtained  were  used  to  construct  on  TRIDAC  a  three-dimensional  model  of 
V.R*725j  complete  vdth  its  control  system,  as  a  replacement  of  the  earlier 
model. 
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After  the  response  and  stability  of  this  model,  which  included  a 
representation  of  the  missile's  autopilot,  had  been  investigated,  the 
simulation  was  extended  to  include  the  guidance  and  homing  system  of 
'V.R.725,  and  this  final  model,  which  required  some  500  amplifiers  and  used 
most  of  the  available  capacity  of  TRIDAC,  was  employed  to  find  the  effect 
of  aerodynamic  characteristics  on  homing  performance.  The  work  on  V.R.725 
forms  the  main  topic  of  this  Note. 

2  MODEL  OF  D.A  TEST  VEHICLE 

2.1  Aerodynamic  characteristics 


A  half-sectional  outline  of  the  D.if  test  vehicle  is  shown  in  Fig,1, 
and  its  principal  dimensions  and  parameters  are  listed  in  Table  1 »  As 
mentioned  in  the  introduction,  information  on  its  aerodynamic  characteris¬ 
tics  was  obtained  from  Ref .2.  These  characteristics  had  been  put  in  mathe¬ 
matical  form  in  Ref. 3,  the  basis  of  the  formulation  being  the  use  of  power 
series  in  the  sines  of  the  partial  incidences  a  and  (3,  defined  in  Fig. 2, 
and  the  standard  S,T),^  representing  the  aileron,  elevator  and  rudder  angles. 
This  scheme  had  been  selected  because  TRIDAC  had  been  provided  with  equip¬ 
ment  designed  to  suit  it,  namoly,  high-grade  hydraulic  seirvo-multipliers 
with  accurate  multiple  potentiometers.  The  proposal  as  envisaged,  however, 
proved  far  too  complicated,  and  was  impracticable  because  the  polynomials 
ropresenting  the  pitching  and  yawing  moments  contained  product  terms  of 
orders  up  to  the  seventh:  these  introduced  awkward  problems  of  scaling 
which  resulted  in  loss  of  accuracy.  A  revised  scheme,  starting  from  the 
assumption  iiiat  the  missile  was  flying  at  constant  height  and  Mach  number 
was  therefore  prepared,  in  which  no  polynomial  was  of  higher  degree  than  the 
third.  The  revised  expression  for  the  pitching  moment  is  quoted  on  Fig«3, 
whore  a  comparison  is  made  with  the  fit  provided  by  the  seventh  order  poly¬ 
nomial.  It  will  bo  seen  that  the  revised  fit  is  certainly  no  worse  than  the 
original  one,  and  in  some  respects  better. 

The  objection  to  both  fits  is  that  they  were  based  on  data  limited  to 
incidences  less  than  22°,  yet  were  employed  to  generate  aerodynamic 
characteristics  at  all  incidences  up  to  30°.  Later  work  showed  the  extrapola¬ 
tion  to  have  been  greatly  in  error,  and  some  of  the  conclusions  roaohod  in 
the  earlier  work  to  be  invalid  for  this  particular  config\iration. 

2.2  Co-ordinate  reference  system 

Fig, 2  indicates  the  system  of  reference  axes  and  notation  employed 

throughout  the  work  as  the  basis  of  the  oquations  describing  the  mathematical 

model,  G,XIZ  are  principal  axes  through  the  mass  centre  G  of  the  missile, 

which  define  its  longitudinal  axis  and  the  planes  through  the  two  pairs  of 

wings.  0  is  the  angle  between  the  velocity  vector  GV  and  the  longitudinal 

axis  and,  therefore,  represents  the  total  incidence.  ^  is  the  angle  between 

tho  XZ  plane  and  the  piano  containing  the  longitudinal  axis  and  the  velocity 

vector.  Hence,  tan  ^  =  v/w,  whore  v  and  w  are  the  sideslip  velocities  in 

the  directions  GY  and  GZ.  The  partial  incidences  a  and  (3  are  used  to  define 

the  quantities  a  =  sin  a  =  w/V  and  p  =  sin  p  =  v/V,  where  V  is  the  speed  of 

tho  missile,  and  the  quantity  5  =  sin  0  is  also  introduced,  so  that 

••  —  — ^  2  — ^ 

a  =  5  cos  p  =  §  sin  (P,  and  0  =  a  +  p  . 

2.3  Equations  of  motion 

The  standard  equations  applicable  to  the  motion  of  a  rigid  body  in 
space,  where  the  motion  is  reforred  to  a  sot  of  moving  axes  which  are 
principal  axes  of  inertia  in  the  body,  were  used,  Tho  equations  are 
listed  in  Appendix  1 .  Since  tho  momonts  of  inertia  about  the  Y  and  Z  axes 
(viz,  B  and  C)  were  assumed  to  be  equal,  the  equations  of  angular  accelera¬ 
tion  were  simplified;  moreover,  the  equation  of  linear  momentum  in  tho 
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direction  of  tho  X  axis  was  replaced  by  a  programmed  variation  of  tho  total 
velocity  V,  or  by  V  =  constant.  What  this  implies  is  indicated  in  Appendix  1 

Tho  quantities  Y, Z,L  M,N,  in  tho  equations  of  motion  are  functions  of 
pressure,  Kach  number  and  Reynolds  number,  as  well  as  a,  3,  etc.  Bosides 
this,  they  depend  to  some  extent  on  p,  q,  r,  and,  except  in  steady  conditions 
on  the  time  derivatives  of  those.  It  was  assumed  however,  that  such  time¬ 
wise  variations  could  be  neglected,  indeed,  in  practically  all  the  work,  tho 
aerodynamic  characteristics  were  assumed  to  be  independent  of  p,  q,  and  r. 

The  wind  tunnel  tests  had  given  no  information  on  rotary  derivatives,  but 
tho  damping  duo  to  them  was  unlikely  to  bo  important  compared  with  the 
artificial  damping  provided  by  tho  autopilot.  Reynolds  number  effects, 
likewise,  wore  unknown  and  had  to  be  neglected,  though  evidence  of  their 
presence  appeared  later  during  flight  trials.  Gravity  forces  wore  neglected, 
since  their  inclusion  in  the  earlier  stages  of  the  work  would  havo  meant 
adding  an  oarth-to-missilo  axis  transformation  to  the  simulation.  After  such 
a  transformation  had  been  added  for  other  reasons  at  a  later  stage,  the 
gravity  terms  in  tho  equations  wore  still  omitted,  since  their  effect  soomod 
likely  to  be  secondary  and  in  tho  nature  of  a  bias  which  might  obscure  the 
main  issue.  Tho  effect  of  gravity  was  included,  however,  in  digital  studies 
undertaken  to  compare  tho  performance  of  tho  model  with  that  of  flight 
rounds.  Since  the  non-linearitios  and  couplings  wore  greater  at  Mach  number 
1*58,  this  value  was  selected  and  adhered  to  in  most  of  tho  work  on  D.4»  An 
altitude  of  40,000  ft  was  chosen  as  suitable,  and  this  fixed  tho  missile 
speed  at  1528  ft/sec. 

2.4  Control  system 

D,4,  like  its  successors  W,1 ,  and  W,2  (V.R.725),  was  fitted  with  an 
autopilot  based  on  rate  feedback  and  designed  to  give  rapid  initial  response 
to  demands  for  rate  of  turn.  Suitable  shaping  networks  wore  employed  to 
modify  the  outputs  of  pitch  and  yaw  angular-rate  gyros  so  that  the  resulting 
lag  in  tho  feedback  path  was  approximately  equal  to  tho  incidence  lag  of  the 
missile.  Since  tho  latter  varied  with  altitude,  however,  provision  vms  made 
to  vary  tho  lag  and  loop  gain  by  switching  at  prodotorminod  altitudes.  The 
shaped  feedback  signals  were  subtracted  from  tho  demanded  rates  of  turn  in 
missile  axes  to  provide  inputs  to  the  control  surface  actuators,  which  were 
driven  by  hydraulic  servo-motors  having  a  natural  frequency  of  about  20  c,p,s 
A  closed  loop  system  also  provided  control  of  roll  rate.  Tho  W.1  and  W,2 
(V.R.725)  autopilots  are  illustrated  in  Pig, 4. 

2.5  Simulator  model 


Pig, 5  is  a  block  diagram  of  the  D,4  simulation.  Tho  main  variables, 

3.  and  3,  which  describe  tho  aerodynamic  characteristics,  were  applied  as 
inputs  to  hydraulic  servo-multipliers  which  formed  thu  polynomial  terms, 

Thoro  wore  insufficient  servos  of  this  typo,  however,  to  generate  tho  control 
surface  contributions  and,  in  consequence,  these  had  to  be  formed  on  electric 
servo-multipliers  of  greatly  inferior  performance.  Por  these  servos  to 
operate  satisfactorily,  it  was  necessary  to  abandon  real-time  working  and 
substitute  a  5:1  time  scale.  Although  this  retarded  progress,  it  was 
advantageous  from  tho  point  of  view  of  accuracy  in  recording  the  results, 
since  tho  full  traverse  of  tho  plotting  tables  could  bo  used.  Digital  check 
solutions  wore  provided  and  tho  simulator  results  shov/ed  good  agreement  vdth 
these . 

2.6  Tests  of  D.4 

As  mentioned  previously,  tho  extrapolation  of  tho  aerodynamic  data  on 
D,4  to  incidoncos  exceeding  22°  was  incorrect;  hence  tho  conclusions  drawn 
from  the  high-incidenco  tests  v/oro  invalid  for  this  configuration.  For  this 
reason,  no  records  of  tho  work  on  D,4  are  included  in  this  Note  other  than 
tho  follovd.ng  general  statement . 
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The  response  of  the  model  to  demands  for  lateral  acceleration  was 
determined  both  in  single  plane,  with  roll  prevented,  and  in  three- 
dimensions,  in  the  presence  of  roll.  Demands  were  applied  in  missile  axes 
and  also  in  axes  which  were  fixed  as  regards  rotation  in  roll  but  which 
coincided  initially  with  the  principal  axes  of  the  missile.  These  were 
termed  "roll-position-stabilised  axes".  Since  the  roll  equilibrium  was 
unstable  when  the  incidence  was  in  either  wing  plane,  but  stable  in  the 
median  planes  (<^  =  45°  or  135°) >  demands  in  roll-position-stabilised  axes 
usually  caused  roll  until  the  wind  vector  lay  in  one  of  the  median  planes. 

As  the  incidence  increased,  instability  developed,  leading  to  divergence  in 
pitch,  or  yaw,  in  a  principal  plane  at  incidences  exceeding  26°  when  roll 
was  prevented,  or  in  combined  pitch  and  yaw  in  a  median  plane,  at  incidences 
exceeding  22°.  This  was  attributable  to  the  markedly  unstable  pitching,  or 
yawing,  moment  which  resulted  from  the  faulty  extrapolation  of  the  data,  a 
fact  which  was  confirmed  by  a  mathematical  analysis  based  on  perturbation 
methods.  When  the  aerodynamic  characteristics  were  artificially  linearised, 
the  instability  disappeared.  In  order  to  prevent  the  instability  from 
affecting  the  rolling  motion,  and  possibly  obscuring  a  genuine  effect, 
linear  characteristics  in  pitch  and  yaw  were  substituted  while  the  non¬ 
linear  rolling  moment  characteristic  was  retained.  There  was  no  evidence  of 
roll  oscillation  under  such  conditions  once  the  missile  had  settled  down 
with  the  incidence  in  a  median  plane,  and  roll  instability  could  be  induced 
only  by  inserting  lags  in  the  roll  control  loop.  These  might  be  thought  of 
as  delays  in  the  generation  of  the  aerodynamic  moment  resisting  roll  -  and 
small  delays  do  occur  owing  to  the  time  taken  for  a  disturbance  of  the  flow 
to  propagate  -  but  the  lags  which  had  to  be  introduced  were  greater  than 
could  be  accounted  for  in  this  way.  An  approximate  analysis  of  the  roll 
loop  was  made,  based  on  the  assumption  that,  since  the  moment  of  inertia  of 
the  missile  in  roll  was  much  less  than  in  pitch  or  yaw,  any  significant 
oscillation  in  roll  must  be  relatively  fast  and,  therefore,  likely  to  have 
little  effect  on  pitch  or  yaw.  The  vrork  on  TRIDAC  had  confirmed  this  theory 
for  the  total  incidence,  6,  had  remained  steady  during  the  roll  oscillations 
Since  an  oscillation  purely  in  roll  implies  p  =  0,  the  conditions  for  oscil¬ 
lation  were  determined  by  analysing  the  roll  loop  on  this  assumption.  The 
results  obtained  completely  confirmed  those  obtained  from  the  simulator, 

3  MODEL  OF  V.R.725 

The  earlier  work  on  the  hybrid  model  based  on  D.4  had  shown  that  more 
extensive  wind  tunnel  tests,  covering  a  greater  range  of  incidence,  were 
essential  if  the  results  were  to  have  any  validity.  In  the  meantime,  the 
D.4,  W.1  and  early  W,2  designs  had  given  way  to  V.R.725  and,  in  the  process, 
the  aerodynamic  configuration,  as  well  as  the  autopilot,  had  been  modified. 

A  completely  now  set  of  wind  tunnel  tests  was  therefore  arranged,  and  these 
were  conducted  in  the  supersonic  wind  tunnel  at  A.R.A. ,  Bedford. 

The  V,R.725  missile  is  a  fixed-wing  cruciform  design  with  four  moving 
control  surfaces  at  the  rear,  A  half-section  on  the  centre  line  is  shown 
in  Eig,1,  and  the  principal  dimensions  and  parameters  of  the  missile  are 
listed  in  Table  1 . 

3.1  Wind  tunnel  tests 


The  ranges  of  incidence  and  Mach  number  covered  in  the  A.R.A.  tests 
were  much  wider  than  those  of  the  early  R.A.E.  tests  on  D.4;  for  example, 
at  Mach  numbers  of  1*2,  1*4,  1*6,  1*8,  2*0,  2*5  and  3*0,  the  maximum 
incidences  reached  in  the  tunnel  v/ere  25°,  25°,  35° >  40°,  40°,  43°  and  43° 
respectively.  The  tests  were  in  two  parts,  the  first  consisting  of  runs 
with  the  control  surfaces  undeflected,  the  second  including  tests  with  one 
control  panel  deflected.  Tests  were  also  made  to  determine  the  interaction 
between  adjacent  control  surface  panels  but,  except  at  the  smaller  Mach 
numbers,  the  interaction  betv/een  panels  was  very  slight.  Refs, 5  to  8 
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discuss  the  organisation  of  these  \Yind  tunnel  tests  and  also  give  the 
results  in  some  detail.  Figs. 6  to  10  summarise  the  aerodynamics  of  V,R,725 
in  the  form  used  in  the  simulation. 

After  the  measured  forces  and  moments  had  been  reduced  to  non-  « 
dimensional  coefficients  by  division  by  a  suitable  factor,  c.g, 

^YpScM^,  the  quotients  were  functions  of  Reynolds  number,  Mach  number,  two 
components  of  incidence  defining  the  position  of  the  wind  vector,  and  three 
angles  representing  the  control  surface  deflections.  As  in  the  caso  of  D,4, 
the  effect  of  Reynolds  number  was  neglected.  It  was  then  necessary  to 
devise  a  scheme  whereby  these  quantities  mi^t  be  generated  on  the  simxilator, 

3.2  Representation  of  aerodynamics 

3.2,1  Tail-fixed  aerodynamics 

Two  methods  v/ere  considered.  The  first  involved_the  expansion  of  the 
force  and  moment  coefficients  as  polynomials  in  a  and  p,  the  method  which 
had  been  used  previously  on  D,4»  TRIDAC  is  well  suited  to  this  method  when 
the  functions  to  be  generated  are  reasonably  linear.  'When  marked  non- 
linearities  are  present,  however,  it  is  much  less  satisfactory.  More  terras 
in  tlie  expansions  are  then  necessary  and,  almost  invariably,  both  positive 
and  negative  signs  occur.  The  associated  constant  multipliers  usually  vaiy 
greatly  in  magnitude  and  present  a  formidable  problem  of  scaling,  with  a 
resulting  loss  of  accuracy. 

The  second  method,  which  was  the  one  actually  employed,  involves  the 
use  of  the  polar  co-ordinates  S  and  <f>  which  have  ali'eady  been  defined.  If 
G-,XY*Z'  is  an  auxiliary  set  of  axes  having  the  X  axis  coincident  with  that 
of  the  missile  but  the  Y  and  Z  axes  rotated  through  the  angle  (p,  so  that  the 
wind  vector  lies  in  the  XZ*  plane,  it  follovjs  from  the  symmetry  of  the 
missile  with  its  control  surfaces  undoflccted  that  the  forces  and  moments 
expressed  in  these  axes  are  periodic  in  40.  The  coefficients  may  therefore 
be  expressed  as  Fourier  sine  or  cosine  series  in  40,  the  coefficients  being 
functions  of  B.  The  general  form  adopted  to  express  the  force  and  moment 
coefficients  in  the  G-,XY*Z'  system,  with  the  control  surfaces  undeflected, 
was  therefore 


TRIDAC  is  well  equipped  to  generate  such  expressions,  since  it  has 
several  highly  accurate  multiple  sine  and  cosine  resolvers  driven  by  hydraulic 
servo  motors.  The  general  term  in  sin  4n0  or  cos  4j^  is  easily  obtained  as 
a  Slim  of  products  of  powers  of  sin  40  and  cos  40.  The  corresponding  coeffic¬ 
ients  in  body  axes  v/ere  obtained  from  the  others  by  resolution  through  the 
angles  0  and  90°  -  0.  The  whole  process  of  generation  and  resolution  was 
accomplished  on  one  servo  by  havirg  hvo  resolver  units  coupled  by  a  gear  box 
of  4j'1  ratio.  The  0  servo  is  described  in  Appendix  2. 

3,2,2  Evaluation  of  &  and  H  functions 

n  n 


The  wind  tunnel  data  were  first  reduced  to  arrays  of  values  of  llie 
force  and  moment  coefficients,  corresponding  to  a  matrix  of  simultaneous 
values  of  6and  0,  itself  the  product  of  serial  values  of  Gand  0,  6  ranging 
in  steps  of  5°  from  zero  to  the  maximum  incidence  attainable  at  each  Mach 
number,  and  0  from  zero  to  45°  in  steps  of  7'^°.  A  greater  range  of  0  was 
unnecessary  by  reason  of  the  symmetry,  althou(^  spot  checks  at  other  angles 
of  roll  v/ere  made. 
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Satisfactory  fits  in  functional  form  were  obtained  in  all  cases  with 

not  more  than  two  terms  in  each  Fourier  series.  The  co-ordinates  of  the 

&  and  H  functions  at  each  e::perimental  value  of  0  were  determined  by 
n  n 

applying  the  method  of  Least  Squares  to  the  appropriate  section  of  the 
data  at  that  incidence.  In  their  final  form,  the  aerodynamic  coefficients 
were  as  follows: 


11 

<D| 

sin  {^  + 

b,  cos 
4 

4) 

c 

m 

=  F^(e) 

+  F2(0)  cos 

14, 

C 

n 

-  ^3(9) 

sin  (1  + 

b^  cos 

4) 

C 

z 

=  F^(9) 

+  Fg(9)  cos 

14 

C 

y 

=  FyCe) 

sin  k4>  {^  + 

b^  cos 

4)  , 

where  the  functions  F  (0)  and  the  coefficients  b  are  introduced  in  place 

of  G  and  H  for  simplicity.  In  these  expressions  b,  ,  b  ,  and  b-  vary 
n  n  M*  J  / 

little  with  5  and  little  error  is  introduced  if  they  are  assumed  constant. 

The  variation  of  the  F  functions  with  incidence  and  Mach  number  is  shown 
in  Pigs ,6  to  9*  The  curves  defining  the  moments  are  particularly  non-linear 
in  0,  and  their  variation  with  Mach  number  is  far  from  regular.  In  fact,  no 
wholly  Satisfactory  method  of  reducing  some  of  these  functions  to  products 
of  functions  of  a  single  variable  was  found.  On  the  simulator,  variation  of 
Mach  number  was  achieved  by  resort  to  linear  interpolation  between  the  values 
of  the  coefficients  at  fixed  Mach  numbers.  Nevertheless,  the  accuracy  of 
representation  of  the  aerodynamics  by  the  method  outlined  was  estimated  to  be 
as  good  as  that  of  the  wind  tunnel  data  from  which  it  was  derived.  Nowhere 
did  errors  exceeding  2  per  cent  of  the  maximum  value  occur  except  where  a 
datum  point  was  completely  out  of  the  pattern,  suggesting  that  the  measurement 
itself  was  wrong, 

3,2,3  Function  generators 

The  F  functions  were  simulated  by  diode-type  function  generators 
v/orking  in  conjunction  with  standard  amplifiers.  These  are  fully  described 
in  Ref ,9,  and  little  further  comment  is  called  for.  Each  generator  was 
initially  set  up  with  a  maximum  error  not  exceeding  1  per  cent  of  full  scale. 
Occasional  checks  were  made  subsequently,  followed  by  adjustaients,  but  the 
maximum  error  did  not  exceed  2  per  cent  of  full  scale  over  periods  of  some 
months.  Seven  function  generators  were  installed  for  the  tail-fixed  aero¬ 
dynamic  coefficients  at  each  Mach  number  and  provision  was  made  for  simula¬ 
tion  at  three  Mach  numbers , 

3«2#4  C ontrol  surf ac e  ae r o dynami cs 

As  might  be  expected,  the  contribution  of  the  tail  surfaces  when 
deflected  consisted  of  a  major  term,  which  was  a  linear  function  of  the 
deflection  of  the  control  surface,  and  a  secondary  part  which  was  a  complica¬ 
ted  function  of  all  the  variables  involved.  Although  the  missile  has  four 
separate  control  panels,  they  are  coupled  together  by  the  autopilot  so  as  to 
provide  only  three  independent  variables,  r\  and  The  transformation 
between  the  two  sets  of  co-ordinates  is  simple. 
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Fig. 10  shows  the  increment  D  on  the  moment  coefficients  caused  by 
deflecting  the  control  panel  cr^  (see  Fig.2) .  The  lateral  aerodynamic  forces 

due  to  the  control  panels  were  assumed  to  be  equal  to  the  moment  divided  by 
a  moment  arm  which  was  taken  as  the  distance  from  the  control  hingo  line  to 
the  centre  of  gravity  of  the  missile.  The  figure  shows  the  data  points,  and 
the  full  lines  represent  the  least-squares  fit  to  these  points  by  a  poly¬ 
nomial  in  a,  p  and  <7^ .  The  effect  of  negative  values  of  <7^ ,  and  of  the 

deflection  of  other  control  panels  may  be  obtained  from  these  curves  by  con¬ 
sideration  of  symmetry,  using  the  relationships  of  Table  2.  Table  2  also 
gives  the  relations  between  the  C7's  and  T)  and  as  defined  in  Fig.2. 

The  poor  control  effectiveness  vdien  (/>  is  near  zero  is  due  to  the  combination 
of  body  incidence  and  control  panel  incidence  viaich  effectively  gives  the 
control  psinel  a  large  incidence.  This  condition  rarely  arises  in  practice. 

It  will  be  noted  that  it  is  accompanied  by  a  relatively  large  and  undesirable 
moment  in  the  other  plane,  giving  rise  to  appreciable  coupling  between  the 
planes.  A  more  detailed  account  of  the  fitting  process,  with  further  results, 
may  be  found  in  Ref. 10. 

The  polynomisils  referred  to  above  were  generated  with  the  aid  of 
electric  servo-multipliers;  not  those  referred  to  in  the  discussion  of  D.2^, 
which  had  teen  discsirded,  but  an  improved  pattern,  described  in  Ref  .11.  The 
effect  of  the  cross-product  terms  arising  from  the  control  panel  deflections 
was  found  to  be  slight;  hence,  in  most  of  the  work  they  were  omitted. 

3.2.5  Damping  derivatives 

The  wind  tunnel  tests  provided  only  the  static  aerodynamic  characteris¬ 
tics  of  the  missile.  Theoretical  estimates  of  the  derivatives  t  ,  m  ,  and  n 

p  q  r 

v/ere  provided  by  Aerodynamics  Department,  R.A.E.,  and  these  were  included  in 
the  model.  The  effect  of  the  latter  two,  in  particular,  was  found  to  be  very 
small  in  comparison  with  the  damping  provided  by  the  autopilot;  thus,  for 
the  most  part,  they  were  omitted  from  the  simulation.  The  effect  of  •ftp,  as 

determined  by  TRIDAC,  is  discussed  in  Section  6. If. 

3 .3  Gravity  forces 

As  in  the  case  of  D.A-,  and  for  reasons  given  earlier  in  this  Note, 
gravity  forces  were  not  included  in  the  simulation,  except  in  a  few  special 
tests. 


3 Variation  of  centre  of  gravity 

Provision  was  made  in  the  digital  programme  for  variation  of  the  posi¬ 
tion  of  the  centre  of  gravity  of  the  missile  resulting  from  the  burning  of 
the  sustainer.  This  causes  a  change  in  the  leverage  of  the  forces  vdiich 
produce  the  moments  about  the  centre  of  gravity.  The  variation  was  assumed 
to  bo  linear  in  time.  Variation  of  the  missile  mass  was  also  included  but 
the  changes  in  the  moments  of  inertia  were  neglected. 

3.5  Variation  of  altitude  and  Mach  number 


To  obtain  the  aerodynamic  forces  and  moments  from  "Hieir  respective 
coefficients,  it  is  necessary  to  multiply  the  latter  by  ^ySpM^  or  -gYScpM^. 
Provision  was  made  to  vary  both  Mach  number  and  altitude,  although  most  of 
the  work  was  done  at  Mach  number  1*6  and  a  simulated  height  of  40,CXX)  ft. 
Better  scaling  was  obtained  by  absorbing  one  M  into  the  F  functions  and 
converting  to  forces  and  moments  by  multiplying  by  KM p/p^^^,  with 

K  =  VY^Pref’  ^^°Pref*  interest  centred  on  large  incidences,  only 
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altitudes  exceeding  that  of  the  tropopause  were  simulated.  The  speed  of 
sound  was,  therefore,  assuimed  constant  and  Mach  niimber  taken  to  be 
proportional  to  the  speed  of  the  missile. 

3 #6  Simulator  model 

A  block  diagram  of  the  arrangement  employed  on  the  simulator  to 
represent  V,R.725  and  its  control  system  is  shown  in  Fig, 11  . 

3*7  Digital  computation 

The  primary  purpose  of  the  digital  computation  was  to  provide  a 
separate  solution  against  which  the  analogue  simulation  could  be  checked. 

It  was,  however,  eventually  used  to  make  comparisons  betv;een  the  mathematical 
model  of  the  missile  and  flights  of  the  actual  missile.  In  the  later  stages 
of  the  programme,  when  it  was  appreciated  that  TRIDAC  was  not  able  to  produce 
miss  distance  values  sufficiently  accurately  for  comparisons  to  be  made 
between  a  missile  with  linear  aerodynamics  and  one  with  non-linear  aero¬ 
dynamics,  the  digital  simulation  was  used  to  produce  the  curves  of  miss 
distance  presented  in  this  Note. 

Because  the  programme  was  originally  intended  to  be  used  for  checking 
purposes  only,  for  which  only  a  few  typical  runs  would  be  required,  some 
computing  speed  was  sacrificed  to  achieve  speed  and  flexibility  of  programm¬ 
ing  and  operation.  Mercury  Autocode  was  used  almost  entirely,  and  all  the 
work  was  done  on  the  Ferranti  Mercury  computer  at  R.A.E.  The  v/ell  tried 
procedure  of  reducing  the  differential  equations  to  a  series  of  first  order 
equations  and  integrating  them  by  a  Runge-Kutta  sub-routine  was  employed. 

The  model  used  was  the  same  as  that  used  on  TRIDAC,  although  the 
TRIDAC  simulation  never  included  all  the  detail  of  the  Mercury  simulation. 

A  series  of  Mercury  programmes  was  written,  but  they  were  all  based  on  the 
first  one  which  was  designed  to  check  the  simulation  of  the  missile  alone. 

They  included  the  provision  of  variable  Mach  number  and  altitude,  a  centre 
of  gravity  position  moving  linearly  with  time,  and  linearly  decaying  mass 
and  inertias.  Slight  modifications  were  made  to  this  programme  to  make  it 
suitable  for  simulation  of  the  flight  rounds.  The  only  major  modification 
necessary  was  when  homing  was  added  to  the  model.  The  provisions  for  varying 
Mach  number,  altitude,  etc,  continuously  were  removed  but  transformations 
between  earth,  missile  and  dish  axis  systems  wore  added.  Radomc  aberation 
and  angular  noise  were  also  included. 

The  F  functions  were  evaluated  by  a  four  point  interpolation  on  a  table 
of  values  of  each  function  at  equally  spaced  intervals  in  6.  The  tables 
were  extended  at  each  end  of  the  range  of  0  so  as  to  avoid  the  time-consuming 
operation  of  making  allowance  for  the  effect  of  being  at  the  end  of  the  range. 

The  transfer  functions  of  the  actuators  v/ere  left  out  in  most  of  the 
work.  The  small  time  constants  would  have  necessitated  a  very  short  inte¬ 
gration  step;  as  it  was,  the  integration  step  was  0*02  sec  for  the  missile 
alone  and  0*01  sec  when  the  homing  head  was  added.  With  the  programme  for 
the  homing  missile  on  which  much  of  the  production  work  was  done,  8  seconds 
of  flight  took  a  little  under  30  minutes  to  compute  -  a  ratio  of  about  210:1. 
This  compares  with  a  time  scale  of  on  TRIDAC.  On  the  analogue  machine 
however,  the  time  taken  to  vary  parameters  and  re-set  initial  conditions  is 
quite  large  compared  with  the  time  of  computation,  whereas  the  opposite  is 
true  of  the  digital  machine.  Hence,  a  more  realistic  figure  for  the  TRIDAC 
computing  ratio  would  be  10:1.  The  time  spent  on  programming  and  eradicating 
faults  must  also  be  reckoned.  This  amounted  to  three  or  four  weeks  in  the 
digital  case,  whereas  the  similar  operation  on  TRIDAC  ran  into  months. 
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4  RESPONSE  OF  V.R.725  MODEL 

4.1  Oeneral 


Many  responses  were  recorded  as  checks  of  the  accuracy  and  symmetry 
of  the  model  and  the  repeatability  of  the  simulator  under  a  wide  variety  of 
conditions.  Most  parameters  of  the  system  v/ere  investigated,  but  Figs.  12 
to  l6  refer  only  to  those  of  most  interest  and.  indicate  the  responses  to 

demands  for  eingular  rates  of  the  missile.  Although  the  actual  demand  in 

each  case  was  a  step  change,  this  was  shaped  into  the  form  indicated  by  a 
simple  lag  before  reaching  the  autopilot,  since  the  latter  is  designed  to 
accept  rate  demands  produced  by  tho  homing  system.  The  effect  of  artificially 

linearising  the  aerodynamics  is  eilso  shown.  For  the  wing-body  contributions, 

this  involved  replacing  the  true  aerodynamic  characteristics  by  straight  lines 
having  the  same  average  slope  over  a  range  of  incidence  from  zero  to  20° 
approximately.  In  the  case  of  the  tail  contributions,  the  major  term  only 
of  the  complete  expression  was  retained  in  each  case.  All  responses  quoted 
are  at  Mach  number  1*6  and  height  40,000  ft,  with  a  fixed  c.g.  position 
corresponding  to  two- thirds  of  the  sustainor  burnt. 

4.2  Responses  in  3in/:;le  plane 

When  the  motion  was  restricted  to  one  principal  plane,  as  with  perfect 
roll  position  stabilisation,  the  responses  were  as  shown  in  Fig, 12  for  tho 
pitch  plane.  The  non-linearity  of  the  pitching,  or  yawing,  moment  caused 
the  double  h\amped  response  in  q,  or  r,  in  contrast  to  the  single  peak  of  the 
response  in  the  case  of  linear  aerodynamics.  The  system  was  completely 
stable  in  both  instances.  The  large  difference  between  the  control  surface 
angles  in  the  two  cases  is  attributable  to  the  choice  of  slope  for  the 
linear  pitching  moment  characteristic.  The  similarity  between  the  plots  of 
incidence,  6,  and  lateral  acceleration,  c,  is  marked,  and  indicates  that  the 
contribution  of  the  control  surfaces  to  the  acceleration  is  small, 

4.3  Responses  in  45°  piano 

Under  equal  demands  for  rate  of  pitch  and  yaw,  tho  missile  remains  in 
the  stable  45°  plane  and  corresponding  parameters  in  the  two  planes  are  equal. 
The  resultant  total  incidence  and  angular  rate  are  in  the  45°  plane  and,  in 
the  case  of  linear  aerodynamics,  are  the  same  as  those  achieved  in  a  roll- 
stabilised  principal  plane  if  the  demand  in  the  latter  case  is  •/2  times  that 
in  both  pitch  and  yaw  in  the  former.  For  this  reason,  responses  in  the  45° 
plane  are  not  illustrated  in  tho  case  of  linear  aerodynamics;  but  Pig, 13 
shows  results  obtained  with  non-linear  wing-body  contributions  which  markedly 
affect  the  operation  of  the  control  system.  In  the  figure  isliich  indicates 
r,^,p,  and  6,  (the  corresponding  curves  in  the  pitch  plane  being  identical) 
equal  demands  r^  and  q^  have  been  selected  to  produce  approximately  the  same 

steady  state  value  of  6  as  in  Fig, 12  for  fully  non-linear  aerodynamics.  The 
demands  are  equivalent  to  a  total  demand  for  0*2  rads/sec,  in  tho  45°  plane, 
but  result  in  much  the  same  steaciy-state  incidence  as  a  single  plane 
demand  foi  0*3  rads/sec.  All  tho  response  curves  are  different,  however, 
and,  in  particular,  it  will  be  observed  that  the  control  surface  angle  in 
Pig, 13,  passes  through  zero  and  has  a  steady  state  value  opposite  in  sign  to 
that  in  Fig,12.  The  explanation  of  this  can  be  seen  in  the  functions  illustra¬ 
ted  in  Pig,11,  which  shows  the  type  of  non-linearity  occurring  in  the  wing- 
body  aerodynamics.  In  the  single  plane  case,  the  angle  kij)  is  eilways  zero  and 
the  functions  F^  and  F^  of  Para. 3. 2, 2  add  to  give  the  total  moment.  In  the 

45°  planes,  however,  40  is  always  l80°  and  the  F^  term  suffers  a  sign  change. 

At  large  incidences,  this  terra  can  be  numerically  larger  than  F^ ,  with  the 

result  that  the  total  moment  becomes  small  and  of  tho  opposite  sign.  The 
necessary  reversal  of  the  control  surface  deflections  is  accomplished  by  the 
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autopilot,  the  missile  rates  of  turn  in  the  steady  state  being  slightly 
greater  than  the  demanded  rates.  In  the  single  plane  case,  the  reverse  is 
true,  the  missile  rate  being  slightly  less  than  that  demanded. 

4-.4-  Response  in  three-dimensions 


The  effect  of  applying  demands  in  roll-position-stabilised  axes  but 
simultaneously  allowing  the  missile  to  roll  is  shown  in  Figs. 14  to  16. 

These  figures  illustrate  the  various  responses  obtained  by  simulating  com¬ 
binations  of  linear  and  non-linear  wing-body  and  tail  aerodynamics  with  the 
same  three  demands  in  each  case.  Here,  the  demand,  denoted  by  where  D 

indicates  "demanded"  and  d  "in  demand  axes",  is  in  the  ZX  plane  of  a  set  of 
demand  axes  having  the  same  origin  and  X  axis  as  those  of  the  missile  but 
not  rolling  with  the  latter.  The  initial  angle  X^,^  between  the  two  ZX 

planes  was  5°  in  order  that  the  missile  should  roll  positively  towards  the 
stable  45°  plane  and  not  be  subject  to  the  uncertainty  existing  at  X^  =  0. 

The  missile's  angular  rates  and  lateral  accelerations  in  its  own  axes  were  . 
resolved  into  the  demand  axes,  giving  the  quantities  q^,  r^,  c^,  and  b^,  p 

being  unaffected  since  no  change  of  roll  axis  occurs. 

The  figures  illustrate  the  motion  which  develops  as  the  missile  rolls 
towards  the  45°  plane,  which  it  eventually  reaches.  The  angular  velocity 
r^  and  the  lateral  acceleration  b^  eventually  decay  to  zero,  whilst  q^  and  c^ 

reach  non-zero  steady  state  values,  that  of  q^  being  close  to  q^^^,  the 

demanded  rate.  The  manoeuvre  is  completely  stable  and  the  double-humped 
feature  of  the  rate  response  under  large  demands  is  again  evident  in 
Figs. 15  and  l6.  The  effect  on  the  response  of  substituting  linear  aero¬ 
dynamics,  wholly  or  in  part,  for  non-lincar  is  shown  in  Figs, 14  and  15  and 
may  be  compared  with  the  result  in  single  plane  (Fig. 12).  Where  only  the 
wing-body  aerodynamic  characteristics  are  non-linear.  Fig, 15  reveals  that 
the  second  hump  of  the  q^  and  r^  responses  is  somewhat  reduced,  as  compared 

with  that  in  Fig,l6,  and  there  is  a  similar  effect  on  the  response  in  roll. 

5  SIMULATION  OF  FLICTiT  CONTROL-ROUNDS 
5,1  Oeneral 


As  mentioned  in  the  Introduction,  a  programme  of  flight  trials  had 
been  envisaged  as  a  possible  addition  to  the  work  on  the  simulator.  For 
various  reasons  which  need  not  be  elaborated,  a  proposal  to  purchase  English 
Electric  control-rounds  for  flight  trials  was  abandoned  but,  during  the 
course  of  the  simulator  exercise,  data  became  available  on  the  performance 
of  similar  rounds  fired  as  part  of  the  development  of  V,R,725.  Some  of  these 
rounds  showed  peculiarities  in  flight  which  merited  further  investigation. 

The  next  sections  of  this  Note  deals  with  this  subject. 

5*2  Flight  rounds 


In  all,  five  of  the  missiles  flown  manoeuvred  satisfactorily  at 
altitudes  and  speeds  relevant  to  the  investigation.  Demands  for  rate  of  turn 
were  in  the  form  of  a  step  through  a  simple  lag.  Various  lags  up  to  2*0  sec 
were  used,  but  the  significant  factor  was  the  nature  of  the  demand  itself, 
that  is,  whether  it  was  a  single  demand  for  acceleration  in  one  wing  plane  or 
equal  demands  in  two  wing  planes  for  manoeuvre  in  a  median  plane.  The  follow¬ 
ing  table  summarises  the  rounds  which  were  of  interest. 
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Round 

033. 

Demand 

10°/3GC 

thro' 

Round 

034. 

Demand 

15°/3ec 

thro  ’ 

Round 

035. 

Demand 

10°/seG 

thro’ 

Round 

036. 

Demand 

15°/3eG 

thro  ’ 

Round 

036. 

Demand 

15°/3ec 

thro' 

Round 

OliJO. 

Demand 

13°/seG 

thro' 

1*5  sec  in  pitch  and  yaw. 


5 #3  Comparison  with  model  results 


The  mathematical  model  was  able  to  reproduce  most  of  the  features  of 
these  rounds  with  a  reasonable  degree  of  accuracy.  Fig.  17  shov/s  a  comparison 
of  the  model  results  with  some  of  the  curves  recorded  by  the  telemetry.  The 
time  interval  shown  covers  the  build  up  of  incidence  in  response  to  the 
demand  quoted  above.  Small  adjustments  to  some  of  the  parameters  in  the 
system  v/ere  needed  for  the  agreement  to  be  as  good  as  that  shown.  ‘The  major 
modification  required  resulted  from  the  fact  that  the  rounds  were  more 
stable  than  the  wind  tunnel  tests  suggested.  A  centre  of  pressure  shift  of 
two  inches  was  required  to  compensate  for  this. 


The  major  disparities  betv/ecn  model  and  missile  occurred  in  the  roll 
response.  Because  there  is  no  position  feedback  in  the  roll  autopilot,  the 
roll  attitude  depends  intimately  on  the  rolling  moment  induced  by  combined 
pitch  and  yaw,  and  the  latter  would  have  to  be  reproduced  accurately  if  the 
roll  response,  as  determined  by  the  simulation,  were  to  agree  with  that  of 
the  rounds.  The  general  character  of  the  roll  response  in  the  simulation 
agreed  well  with  that  observed  in  the  rounds  except  for  two  details  which  may 
well  be  related,  but  for  which  no  explanation  has  been  found.  Fig.l8  shows 
extracts  from  the  records  of  roll  rate  gyro  output  of  four  of  the  rounds,  and 
also  gives  some  indication  of  the  range  of  speed  and  altitude  covered.  In 
the  rounds  where  the  demand  was  for  rate  of  turn  in  one  wing  plane  only,  the 
induced  rolling  moment  made  the  missile  roll  in  an  attempt  to  reach  the  stable 
state  of  equal  incidences  in  two  planes.  In  the  transient  so  formed  there 
were  small  pealcs  in  the  records  from  the  rounds,  which  were  not  reproduced  by 
the  simulation.  This  is  shown  for  Round  03A  in  Fig.l8.  The  other  three 
rounds,  records  of  which  are  shown  in  this  figure,  are  those  in  which  the 
stable  state  with  equal  incidences  in  two  wing  planes  should  have  been  attained 
without  any  disturbances  in  roll;  but,  as  can  be  seen,  each  one  developed  an 
oscillation.  A  considerable  amount  of  time  was  spent  in  trying  to  make  the 
model  oscillate  in  a  similar  manner,  so  that  the  cause  of  this  oscillation 
might  be  understood,  but  no  real  success  v/as  achieved.  The  telemetry  records 
of  Round  035  were  particularly  good,  showing  the  oscillation  to  be  of  about 
3  c/s  frequency  with  a  modulating  frequency  of  about  0*3  c/s.  A  detailed 
analysis  of  the  records  at  various  points  in  the  roll  loop  was  made  by  the 
manufacturers.  Sections  of  the  model  were  checked  independently  and  all  but 
the  aerodynamic  rolling  moment  appeared  correct.  The  work  on  TRIDAC  is 
described  in  the  next  section. 

6  INVE3TI&ATI0N  OF  ROLL  STABILITY 

6.1  General 


¥/ith  the  aerodynamic  and  control  system  data  appropriate  to  M  =  1*6, 

H  =  40,000  ft,  which  were  adopted  as  representative  values,  there  was  no 
evidence  on  the  simulator  of  roll  oscillation  at  any  incidence  up  to  30^. 
Attempts  were  made  to  induce  oscillations  by  various  methods,  viz:  (l)  intro¬ 
duction  of  a  small  lag  in  the  forv/ard  path  of  the  roll  loop,  i.e.,  in  L/A, 

(2)  introduction  of  a  small  lag  in  the  feedback  path  of  the  roll  loop,  i.e., 
in  together  with  a  small  lag  in  thedemandfor  roll,  i.e.,  in  (3) 

introduction  of  backlash  and  dead-spacG  in  the  control  surface  displacements. 
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(4)  changes  in  gains  and  time  constants  in  the  autopilot  loops  of  up  to 
±10  per  cent,  (5)  simulation  of  non-linear  control  surface  aerodynamics. 

The  results  of  these  experiments  are  described  below. 

6.2  Investigation  of  roll  loop 

Since  the  circular  frequency  of  the  roll  oscillation  observed  in 
flight  had  been  recorded  as  20  rads  per  second,  approximately,  the  outer 
loop,  though  the  couplings  with  pitch  and  yaw,  was  broken  at  the  demand 
for  roll,  i.e.  at  the  L^<}>  input,  and  amplitude  and  phase  measurements  were 

made  at  this  frequency.  The  values  obtained  were  0*59  and  -305*5°  respec¬ 
tively,  In  order  to  produce  instabilily  in  roll  due  to  the  outer  loop, 
therefore,  it  would  be  necessary  to  introduce  a  gain  of  about  1*7  and  a 
phase  lag  of  54*5°  into  the  outer  loop.  Assuming  that  the  value  of  was 

correct,  it  would  be  possible,  on  a  linear  basis,  to  do  this  by  the  intro¬ 
duction  of  a  second  order  lag  having  a  damping  ratio  of  0*5  and  a  time 
constant  of  41  m,s.,  or  by  a  phase  delay,  1-sT/l+sT,  with  a  time  constant  T 
of  26  m,s.,  together  with  a  gain  of  1*7»  For  a  simple  lag  of  70  m.s.,  the 
gain  Inquired  would  be  about  2*9.  With  a  symmetrical  missile,  none  of  these 
possibilities  seemed  feasible.  It  was  assumed,  therefore,  that,  apart  from 
a  possible  variation  in  the  magnitude  of  L,,  any  lag  occurring  in  the  system 
must  be  in  the  roll  loop  itself.  ^ 

6 .3  Effect  of  lag  in  L/A 

In  order  to  investigate  the  effect  of  lags  within  the  inner  roll  loop, 
equal  demands  for  q  and  r  were  provided  in  roll  position  stabilised  axes  so 
that  the  missile  took  up  a  position  in  the  45°  plane.  The  demands  were  such 
that  the  steady  state  value  of  5  was  0*462,  corresponding  to  an  incidence  of 
27*5°,  A  lag  of  36  m.s.  in  L/A  then  produced  a  continuous  oscillation  in  p, 
and  increasing  the  time  constant  of  the  lag  caused  instabilily  in  roll.  The 
frequency  of  the  oscillation  was  3*4  c.p.s.,  which  was  of  the  right  order. 

The  effect  of  changing  tiie  value  of  is  shown  in  Fig. 19.  Similar  results 

were  obtained  at  altitudes  of  25,000  ft  and  53,500  ft,  the  lag  required  being 
reduced  at  the  lower  altitude.  Further  results  were  obtained  for  Mach 
numbers  of  1*4  and  2*0,  and  these  showed  that  the  lag  required  was  slightly 
less  at  the  higher  Mach  number.  At  the  normal  value  of  L.,  the  smallest  lag 
required  to  cause  oscillation  was  25  m.s.,  at  Mach  number^2*0,  and  altitude 
25,000  ft.  The  frequency  of  oscillation  was  5*5  c.p.s.  With  a  gain  of  2*0 
on  L^,  this  lag  was  reduced  to  12  m.s.  On  the  assumption  that  the  value  of 

L.  used  on  the  simulator  was  correct,  the  lag  required  to  produce  roll 

oscillation  was  considered  to  be  more  than  was  likely  to  arise  in  the  roll 
autopilot  loop. 

The  results  obtained  in  this  section  of  the  work  were  checked  by 
perturbation  theory,  good  agreement  being  obtained.  The  digital  check  solu¬ 
tions  gave  similar  results, 

6.4  Effect  of  roll  damping  derivative  6 


As  was  expected,  the  introduction  of  a  value  of  (which  had  been 

estimated  to  be  -6*3)  had  a  stabilising  effect.  V/ith  a  positive  •6^  of  the 

same  magnitude,  the  stabilily  was  reduced,  but  as  much  as  16  m.s.  lag  was 
still  required  to  produce  continuous  oscillation  at  M  =  1*6'  and  H  =  40,000  ft. 
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6.5  Effect  of  lag  in 


It  is  unlikely  that  any  lag  in  could  occur  unless  it  had  an  aero¬ 
dynamic  origin,  but  lags  in  are  possible  if  imperfections  exist  in  the 

control  system  components  and  in  the  roll-rate  gyro.  It  was  not  possible, 
however,  to  start  an  oscillation  in  roll  by  means  of  such  lags  alone.  With 
a  lag  of  50  m.s.  in  for  example,  a  lag  of  10  m.s.  in  was  necessary 

to  steurt  oscillation  at  an  incidence  of  30®.  Various  combinations  of  lags 
in  and  L^(f>  were  explored,  the  i:osults  being  shown  in  Fig, 20, 

6 ,6  Theoretical  estimate  of  aerodynamic  lags 


It  had  been  suggested  that  lags  in  the  control  system  might  occur  owing 
to  the  finite  time  required  for  aerodynamic  loads  to  reach  steady  values  after 
a  disturbance.  Tobak"^,  who  used  indicial  functions  to  study  unsteady  airflow, 
concluded  that  the  time,  t,  required  for  loads  on  lifting  surfaces  to  atteiin 
steady  values  after  a  sudden  change  of  flow  is  equal  to  the  chord  of  the  lift¬ 
ing  surface  divided  by  the  velocity  at  which  the  distxirbance  passes  downstream, 
i.e.. 


t 


L 

u-a 


where  t 


time  in  seconds 


L  =  streamwise  chord  of  lifting  surface  (i.e.  wing  chord)  or, 

eO-ternatively,  distance  from  wing  leading  edge  to  tail  leading 
edge 

u  =  free  stream  velocity 
a  =  speed  of  sound  in  free  stream 


In  the  case  of  Round  035>  the  Mach  number  and  altitude  when  the  roll 
oscillations  commenced  indicated  a  lag  of  about  17  m.s.  for  the  wing- tail 
combination,  and  a  lag  of  8*5  m.s.  for  the  wing  alone.  These  times  are  too 
short  to  account  for  the  oscillation,  as  the  work  described  above  shows, 

6.7  Effect  of  control  equipment  non-linearities 

To  determine  the  contribution,  if  any,  of  possible  non-lino arities  in 
the  control  geeir  to  roll  instability,  dead-space  and  backlash  were  simulated, 
the  former  to  represent  overlap  in  the  hydraulic  control  valves  and  the 
latter  to  represent  mechanical  imperfections  in  the  control  surface  drives. 

Dead-space 

It  was  found  necessary  to  simulate  control  valves  in  which  the  overlap 
was  25  per  cent  of  the  total  travel  before  an  oscillation  in  roll  of  tho 
right  amplitude  was  produced.  In  practice,  such  overlaps  are  nominally  zero 
and,  in  any  case,  would  bo  very  small;  hence  dead-space  was  ruled  out  as  a 
possible  cause  of  instability. 

Backlash 

A  realistic  amount  of  backlash,  viz;  0*1°  in  a  total  travel  of  30®  was 
simulated  but  did  not  cause  any  noticeable  decrease  in  roll  stability.  An 
increase  of  backlash  beyond  this  amount  could  not  be  justified;  hence,  it 
was  concluded  that  a  non-linearity  of  this  type  could  be  discounted. 
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6.8  Effect  of  loop  gains  and  timo  constants 

For  variations  in  the  gedns  and  time  constants  in  the  autopilot  loops 
between  ±10  per  cent,  no  noticeable  effect  on  roll  stability  was  observed. 

6.9  Effect  of  control  surface  aerodynamic  non-linearities 

The  non-linearities  in  the  control  surface  aerodynamic  characteristics, 
when  simulated,  had  a  stabilising  effect  on  the  roll  loop  for  incidences  up 
to  30®. 

6.10  Effect  of  missile  asymmetries 

All  the  tests  mentioned  so  far  were  made  on  the  assumption  that  the 
missile  was  completely  symmetrical  in  both  the  pitch  and  yaw  pianos.  To  dis¬ 
turb  the  symmetry,  various  biases  were  added  to  the  force  and  moment  equations 
but,  for  what  were  considered  to  be  reasonable  values,  no  noticeable  effect  on 
roll  stability  was  produced. 

Also,  to  upset  the  symmetry  of  the  system,  gravity  components  were 
included.  Although  these  appeared  to  decrease  the  stability,  and  an  oscilla¬ 
tion  for  a  limited  time  was  observed,  its  amplitude  was  very  small  compared 
with  that  of  the  flight  rounds. 

6.11  Conclusions  on  roll  instability 

V/ith  the  V.R.725  missile  and  control  system  simulated  as  described,  no 
significant  roll  oscillation  could  be  produced  by  the  Introduction  of  what 
were  considered  to  be  reasonable  lags  and  changes  in  gains,  time  constants 
and  other  parameters.  It  was  concluded,  therefore,  that,  if  the  oscillations 
in  roll  observed  in  flight  represented  a  genuine  instability,  its  cause  must 
be  aerodynamic  -  probably  a  transient,  or  unsteady,  phenomenon  which  had  not 
been  observed  in  the  wind  tunnel  and  which,  therefore,  had  not  been  simulated. 

7  I^WESTI&ATION  OF  HOMINC  PERFORMANCE 


7.1  System  adopted  for  study 

The  ultimate  criterion  of  perfomance  in  a  missile  aimed  at  a  target 
is  lethality.  If  fuse  and  warhead  performances  are  not  to  be  considered, 
miss  distance  is  generally  accepted  as  the  criterion  of  performance.  It  was 
a  natural  extension  of  the  work, therefore,  to  seek  to  determine  the  effects 
of  the  aerodynamic  non-linearities  and  cross-couplings  on  the  miss  distance 
achieved  by  the  missile.  The  simulation  was  extended,  therefore,  to  include 
the  V.R.725  homing  head  with  its  servo  dish  and  feedback  loops. 

7*2  Details  of  extension  to  simulation 


Two  axis  transformations  had  to  bo  constructed.  The  first  of  these  v/as 
a  three-dimensional,  three  gimbal  system  to  transform  the  assumed  target 
velocities  from  earth  to  missile  axes.  The  second  was  a  two  gimbal  system 
to  transfonn  the  relative  velocities  of  missile  and  target,  and  the  missile 
rates  of  turn,  from  missile  to  dish  axes.  These  transformations  were 
accomplished  with  the  aid  of  five  hydraulic  servo-resolvers  which  are  part 
of  the  equipment  of  TRIDAC  provided  specially  for  such  purposes.  The  servos 
may  be  used  either  as  velocity  servos  with  tachometer  feedback,  or  as  position 
servos  with  position  potentiometer  feedback.  At  the  start  of  the  work,  the 
former  arrangement  was  employed,  but  it  was  abandoned  later  in  favour  of  the 
latter,  on  account  of  the  errors  caused  by  a  small  amount  of  drift  in  the 
servos.  In  the  early  stages  also,  the  relative  velocities  of  missile  and 
target  were  integrated  in  missilo  axes  before  transformation  into  dish  axes. 
This  resulted  in  the  appearance  of  considerable  noise  superposed  on  the 
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components  of  relative  position  in  dish  axes,  owing  to  the  limited.,  resolution 
of  the  potentiometers .  The  trouble  was  overeorae  by  integrating  the  relative 
veloeities  in  dish  axes  after  their  transformation  from  missile  axes. 

Integration  in  dish  axes  brought  the  further  advantage  that  tho  sealing 
could  be  improved.  Usually,  in  a  homing  manoeuvre,  the  dish  points  very 
nearly  at  the  target  throughout,  so  that  the  lateral  eomponents  of  target 
position  in  dish  axes  are  always  small.  Only  the  axial  eomponent  varies 
greatly  in  magnitude.  Measured  in  missile  or  earth  axes,  however,  all 
three  eomponents  of  relative  position  are,  in  general,  large  Tirtien  the  missile 
and  target  are  far  apart,  and  small  when  they  are  near.  Consequently,  the 
voltages  representing  the  lateral  eomponents  of  miss  distanee  in  missile  or 
earth  axes  are  small  in  relation  to  their  maximum  values  when  tho  missile 
and  target  aro  near,  whereas  this  is  not  true  of  tho  eomponents  in  dish  axes. 

In  order  to  form  error  signals  representing  the  radar  output,  tho  tv/o 
lateral  eomponents  of  target  position  in  dish  axes  were  divided  by  tho  third 
component  measured  along  the  dish  axis  which  is  virtually  the  range.  The 
division  proeess  was  aecomplishod  with  tho  aid  of  an  eleetric  servo- 
multiplier.  •  The  V.R. 725  homing  system  uses  two  sueh  eiror  signaDs  as  inputs 
to  the  dish  servo  drives,  and  also  as  the  demands  for  missile  rate  of  turn 
to  the  pitch  and  yaw  autopilots.  The  simulator  model,  when  extended  to 
include  the  homing  systean  and  target  motion  as  indicated  above,  is  shown 
diagrammatieally  in  Pig. 21,  with  the  detailed  representation  of  the  homing 
head  in  Pig. 22. 

7.2.1  Meastirement  of  miss  distanee 

As  was  stated  earlier,  the  assumed  eriterion  of  homing  performanee  was 
the  miss  distanee.  Various  methods  of  measuring  miss  distance  were  considered 
and  tried.  The  objection  to  foming  the  sum  of  the  squares  of  the  displaoe- 
ments  and  extracting  the  square  root,  i.e.,  the  range,  is  that  the  magnitude 
of  the  latter  varies  so  much  that  good  scaling  is  impossible.  As  an  alterna¬ 
tive,  one  may  form  the  predieted  miss  distanee,  i.e.,  the  miss  that  would 
result  if,  subsequent  to  the  instant  of  predietion,  tho  relative  motion  were 
linear  at  eonstant  voloeity.  This  quantity  has  the  advantage  that  it  varies 
more  slowly  than  the  range  but,  in  three  dimensions,  its  formation  is  sotio- 
v/hat  complieated  and  the  proeess  is  liable  to  error  unless  highly  aecurate 
equipment  is  available.  In  the  end,  it  was  found  best,  when  studying 
approaehing  targets  and  also  for  eomparative  purposes,  to  fom  the  resultant 
of  the  two  lateral  eomponents  of  range  in  missile  axes  and  reeord  it  over  the 
last  few  seeonds  at  the  end  of  the  engagement  in  order  to  determine  the 
minimum.  A  further  convenient  assumption  made  when  miss  distanoes  were  being 
recorded  for  eomparativo  purposes  was  that  the  minimum  range  occurred  when 
its  component  along  the  missile  longitudinal  axis,  was  zero.  This  is  not 
strietly  aeeurate,  but  it  is  a  good  approximation  when  dealing  with  nearly 
head-on  targets.  The  miss  distance  was  then  recorded  by  feeding  the  voltage 
representing  the  root-raean-square  of  the  two  lateral  components  of  range  in 
missile  axes  to  the  initial-condition  input  of  an  integrator,  the  signal 
input  to  liie  integrator  being  grounded.  When  the  longitudinal  eomponent  of 
range  went  through  zero,  the  integrator  was  started  auto matieally,  so  that, 
its  output  was  held  at  the  miss  distance  as  defined  above.  This  output  was 
then  displayed  on  a  digital  voltmeter.  Besides  miss  distanee,  the  other 
variables  of  chief  interest  wore  reeorded  eontinuously  on  a  conventional 
servo-driven  plotting  table,  and  on  a  12-ehannel  ultra-violet  roeorder.  The 
latter  had  the  advantage  of  good  response,  whereas  the  former  eould  be  used 
for  eomparing  the  results  of  two  engagements  by  direet  superposition  of 
reeords.  The  target  was  always  assumed  to  be  travelling  at  1000  ft/see  and 
the  missile  at  a  Maeh  number  of  1*6  at  40,000  ft  altitudo. 

At  a  range  of  300  ft,  the  homing  system  was  assumed  to  be  saturated 
by  the  signal,  and  a  diodo  switeh  was  ineluded  to  stop  the  range  servo  \:Aien 
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the  longitudinal  component  of  range,  measured  in  dish  axes,  reached  this 
value.  As  an  alternative,  the  switch  was  arranged  to  earth  the  inputs  to 
the  dish  actuators  and  autopilot,  thus  reducing  the  demand  to  zero.  The 
difference  in  recorded  miss  distance  was  usually  negligible,  so  the  precise 
assumption  made  did  not  appear  to  matter.  Miss  distance  was  also  found  to 
be  insensitive  to  a  variation  in  the  cut-off  point  between  500  and  500  ft . 

7,2,2  Acceleration  limiter 

Although  interest  in  the  work  was  primarily  in  the  effects  of  high 
incidence,  the  aerodynamic  data  did  not  extend  beyond  35° •  A  lateral 
acceleration  limiter  was  installed  to  contain  the  incidence  within  any 
desired  range.  It  took  the  form  of  a  diode  switch  which,  when  the  accelera¬ 
tion  reached  the  limit,  operated  to  reduce  the  inputs  to  the  pitch  and  yaw 
autopilots  to  one  hundredth  of  their  normal  values.  This  limit  was  normally 
set  to  operate  when  the  incidence  reached  about  30°.  Whenever  possible, 
homing  runs  in  which  this  limit  was  not  reached  were  favoured,  as  the  effect 
of  the  limiter  confused  the  results. 

It  should  be  mentioned  that,  throughout  the  simulation,  for  reasons 
already  given,  the  time  scale  employed  was  five  times  real  time, 

7,3  Targets  investigated 

Homing  was  studied  against  fixed  targets,  straight-flying  targets,  and 
targets  making  manoeuvres  such  as  weaves  and  steady  turns.  In  performance 
studios,  the  weaving  target  is  often  considered  to  be  the  most  difficult  one 
to  intercept.  Large  miss  distances  are  obtained  because  the  navigation 
system  is  not  well  suited  to  the  requirements  in  this  case.  Weaving  target 
manoeuvre  was  not  considered  very  appropriate  in  the  present  investigation 
however,  because,  although  it  was  tried,  and  resulted  in  large  miss  distances 
on  occasion,  the  missile  usually  did  not  attain  large  incidences  during  its 
homing  run.  Since  interest  centred  on  effects  occurring  mainly  at  large 
incidences,  a  steadily  turning  target  v/as  more  suitable. 

It  is  only  in  the  final  stage  of  an  engagement  that  incidences  are 
large,  so  this  was  the  only  part  of  the  homing  attack  considered.  The 
standard  engagement  was  one  in  which  the  missile  and  target  were  assumed  to 
be  approaching  at  constant  speed  at  an  altitude  of  4-0,000  ft  until  they  were 
within  20,000  ft.  The  target  was  then  assumed  to  begin  a  steady  2g  turn, 
this  was  the  starting  point  of  the  simulation.  Target  motions  were  produced 
by  resolving  the  total  velocity,  V^,  of  the  target,  assumed  to  be  a  point, 

into  components  u^^^  and  or  u^^  and  in  fixed  earth  axes. 

8  EFFECT  OF  AERODYNAMICS  ON  MISS  DISTANCE 
8,1  General 


In  the  earlier  part  of  the  programme,  the  method  of  assessing  the 
effect  of  aerodynamic  non-linearities  and  cross-couplings  had  been  to  compare 
the  performance  of  the  missile  with  that  of  a  hypothetical  missile  identical 
with  it  except  for  the  absence  of  non-lineardties  and  cross- couplings  from 
its  aerodynamics.  In  fact,  however,  when  it  came  to  comparing  miss  distances, 
it  was  found  that  TRIDAC  did  not  generate  the  latter  accurately  enough  for  the 
comparison  to  bo  meaningful.  To  obtain  the  required  accuracy,  it  was  necessary 
to  resort  to  digital  simulation,  using  the  MERCURY  autocode  programme,  which 
had  been  written  primarily  to  check  TRIDAC.  As  each  run  took  about  half  an 
hour  on  MERCURY,  the  range  of  conditions  investigated  had  to  be  curtailed  and, 
thus,  a  number  of  parameters  which  might  have  been  varied  had  to  be  fixed. 
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The  digital  results,  presented  in  Pigs ,23  to  27,  are  designed  to  show 
the  effects  of  the  various  forms  of  non-linearity  and  cross-coupling,  as 
well  as  the  interactions  between  them.  Non-linear  characteristics  were 
replaced  by  straight  lines  to  show  the  effect  of  the  non-linearitieB,  and 
the  rolling  moment  due  to  combined  pitch  and  yaw  was  suppressed  to  indicate 
its  effect. 

8.2  Homing  in  single  plane 

The  initial  range,  i.e.,  the  range  at  which  the  target  started  its 
manoeuvre,  was  varied  and  miss  distanoe  was  measured  \^4lilst  the  missile  was 
constrained  to  fiy  with  incidence  in  one  wing  plane  only.  The  effect  of  non- 
linearities  in  pitch  (or  yaw)  is  shown  in  Fig. 23.  It  will  be  noted  that  the 
miss  distance  is  always  greater  when  non-linearities  are  present.  As  the 
initial  range  is  decreased,  the  missile  is  forced  to  manoeuvre  at  progres¬ 
sively  higher  incidences  to  achieve  the  necessarily  higher  accelerations 
(Pig,2A)*  The  disparity  between  the  linear  and  non-linear  results  increases 
to  the  point  where  the  acceleration  limiter  operates.  At  still  shorter 
initial  ranges  the  missile  has  little  time  in  which  to  manoeuvre,  so  that 
the  miss  distance  is  largely  determined  by  the  lateral  displacement  of  the 
target  during  the  engagement.  The  change  of  sign  of  miss  distance  indicates 
that  beyond  a  certain  range  the  missile  over-anticipates  the  target's  motion 
and  passes  the  target  on  the  opposite  side. 

8.3  Homing  in  three  dimensions 

The  effects  of  the  in-plane  non-linearities  may  now  be  compared  with 
those  arising  from  non-linearities  and  cross-couplings  in  three  dimensions. 
Fig. 25  compares  results  at  two  ranges,  viz,  1^,000  ft  and  20,000  ft,  where 
the  peak  incidence  reached  during  the  engagement  ranges  from  21^  to  25^  and 
from  17°  to  18°,  respectively.  A  new  parameter  introduced  by  the  third 
dimension  is  the  roll  attitude  of  the  missile  relative  to  the  plane  in  which 
the  target  is  manoeuvring.  Since  the  roll  autopilot  serves  only  to  damp  the 
rolling  motion,  the  missile  is  free  to  assume  any  roll  position  during  a  run. 
Hence,  the  parameter  chosen  to  represent  this  new  degree  of  freedom  was  the 
initial  roll  angle  ihe  definition  of  which  can  most  easily  be  understood 

by  reference  to  Pig. 28.  It  will  be  recalled  that  the  rolling  moment  due  to 
combined  pitch  and  yaw,  which  alone  causes  the  rolling,  is  zero  for  (f>  =  0, 

45,  90,  135,  etc,  degrees,  these  angles  representing,  alternately,  unstable 
and  stable  roll  positions.  During  the  early  part  of  the  engagement,  before 
the  missile  has  rolled,  the  angles  Xp  and  (f>  will  be  nearly  equal,  (not  quite 

U- 

equal,  because  the  responses  in  the  two  planes  differ  slightly  since  the  dish 
is  mounted  on  two  gimbals,  one  of  which  is  inside  the  other);  but,  as  its 
incidence  increases,  the  missile  will  tend  to  roll  so  as  to  make  <f>  approach 
one  of  its  stable  values.  Up  to  a  point,  the  nearer  the  initial  position  to 
one  of  unstable  equilibrium  in  roll,  {<f>  =  0,  90,  180,  ...  degrees),  the  more 
violent  is  the  rolling  motion  and,  as  will  be  seen  from  Pig. 25,  there  is  a 
corresponding  increase  in  miss  distance.  Fig. 25  also  indicates  that  the  non- 
linearities  of  the  pitch  and  yaw  aerodynamics  have  their  effect;  for  the 
increase  in  miss  distance  does  not  occur  if  either  the  rolling  is  suppressed 
or  the  pitch  and  yaw  aerodynamic  characteristics  are  linearised.  The 
asymmetcy  in  the  curves  is  due  to  the  asymmetry  in  the  dish  gimbal  system. 

Although  to  suppress  rolling  entirely  is  somewhat  artificial,  it  may  be 
regarded  as  changing  the  control  system  to  one  of  perfect  roll  position 
stabilisation.  The  results  suggest,  therefore,  that  roll  position  stabilisa¬ 
tion  is  preferable  to  stabilisation  of  roll  rate,  at  any  rate  for  the  type 
of  homing  manoeuvre  investigated. 
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9  OTHER  EFFECTS 

9.1  General 


It  ia  important  to  know  whether  the  large  peaks  in  miss  distance 
referred  to  in  the  previous  sections  are  reduced  by  some  of  the  phenomena 
present  in  the  missile  system,  which  have,  so  far,  not  been  mentioned.  The 
effects  of  noise  and  of  a  steady  non-zero  rate  demand  were  investigated. 

9.2  Effect  of  angular  noise 

In  order  to  assess  the  effect  noise  would  have  on  the  curve  of  miss 
distance  vsiriation  with  initial  roll  angle,  several  runs  were  made  with 
angular  noise  included.  Noise  which  was  isiiite  within  the  bandwidth  of  the 
homing  head  and  of  spectral  density  lO"?  rad^/rad/sec  was  added  to  the 
error  inputs  to  the  dish  system.  The  two  sources  of  noise  were  uncorrelated. 
Pig, 27  shows  the  miss  distances  obtained  with  an  initial  range  of  14,000  ft, 
and  compares  them  with  the  noise-free  results.  It  can  be  seen  from  the  Pig. 
that  the  effect  is  not  one  of  smoothing  out  the  curve  but  of  increasing  the 
miss  distance  overall. 

9.3  Effect  of  steady  roll 


Pig, 26  compares  the  miss  distances  obtained  under  a  steady  roll-rate 
demand  of  0*25  rad/sec  with  those  for  zero  roll  rate  demand  already  shown  in 
Pig, 25.  The  initial  range  in  this  case  was  20,000  ft.  Again,  the  effect  is 
not  to  smooth  out  the  curve  but  rather  to  cause  a  general  increase  in  miss 
distance.  It  is  understood  that  V.R.725  has  a  steady  demanded  rate  of  roll 
of  this  order  to  reduce  the  effects  of  possible  biases. 

10  SIMULATOR  ACCURACY 

10.1  General 

The  simulation  entailed  a  great  deal  of  duplication  in  that  identical 
pitch  and  yaw  planes  were  represented,  gravity  being  omitted.  Hence,  a 
valuable  check  on  correct  operation  could  be  made  by  comparing  responses  in 
separate  planes  to  equivalent  demands;  further,  for  either  plane,  checks 
of  symmetry  of  response  for  symmetry  of  demand  were  possible  on  all  para¬ 
meters.  Thus,  if  four  responses  of  identical  shape  could  be  achieved, 
together  with  good  repeatability,  oonfidence  that  the  simulator  was  function¬ 
ing  correotly  would  be  increased.  If  these  responses  then  matched  those 
obtained  digitally  to  a  satisfactory  degree,  confidence  that  the  system 
simulated  was  correct  would  follow.  Similar  comparisons  could  also  be  made 
in  the  45°  planes  to  provide  additional  checks . 

10.2  Missile  without  homing  system 

Excellent  repeatabllily  of  all  parameters  associated  with  response  to 
missile  rate  demands  was  obtained,  together  with  good  symmetry.  As  shown  by 
Pig, 29,  however,  agreement  with  the  digital  check  solution  was  less  satis¬ 
factory  in  some  places,  the  greatest  error  recorded  being  about  3  per  cent 
of  the  full  scale  vsilue.  As  far  as  the  investigation  of  stability  was  con¬ 
cerned,  these  errors  could  be  Ignored,  although  considerable  effort  was  made 
to  reduce  them,  without  success.  Their  presence  was  a  reminder  that  TRIDAC 
was  not  built  for  absolute  measurements  on  complex  systems,  and  presaged  an 
increasing  divergence  from  the  digital  solutions  as  complexity  grew. 

10.3  Homing  missile 

The  homing  system  was  simulated  in  considerable  detail,  and  a  great 
increase  in  the  amount  and  complexity  of  the  equipment  resulted.  Where 
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possible,  individual  components  of  the  homing  simulation  were  checked  with 
the  aid  of  calculated  responses,  and  tests  for  symmetry  were  made  similar  to 
those  applied  to  the  missile.  An  acceptable  standard  was  obtained  in  all 
oases,  after  which  the  system  was  coupled  up. 

A  deterioration  in  the  accuracy  of  the  simulator  was  immediately 
apparent  in  homing  manoeuvres.  After  a  fairly  lengthy  process  of  elimination 
of  unreliable  components,  short  term  repeatability  of  responses  was  achieved, 
longer-term  random  variations  remaining.  These,  however,  were  less  than  the 
general  asymmetry  of  the  responses  to  symmetrical  inputs  (applied  in  single 
planes  or  with  X_  =  45°) •  It  was  significant  that  the  accuracy  with  which 

VT 

the  various  response  curves  matohed  the  digital  solutions  had  decreased. 

Methods  used  to  try  to  improve  the  accuracy  of  the  system  and  to  pin¬ 
point  the  chief  factors  leading  to  an  accumulation  of  errors  are  described 
in  Appendix  3*  Ultimately,  the  errors  in  response  were  accepted  and  emphasis 
was  placed  on  miss  distance  measurement.  As  expected,  long  term  repeatability 
and  symmetry  in  the  measurements  were  poor,  but  it  was  found  that  the  average 
of  four  values  of  miss  distance,  obtained  under  symmetrical  conditions  of 
target  manoeuvre,  often  agreed  well  with  the  digital  figure,  see  Figs .30  and 
31.  It  was  considered,  therefore,  that  useful  information  could  be  obtained 
from  the  simulator  on  the  general  effect  of  system  changes  on  miss  distance. 

11  CONCLUSIONS 


Marked  non-linearities  and  cross- couplings  occur  in  the  aerodynamic 
characteristics  of  cruciform  missiles  of  conventional  design  -  of  which 
V.R.725  may  be  regarded  as  typical. 

The  most  important  characteristic  is  a  rolling  moment  which  increases 
with  incidence  and  which,  if  the  missile  is  not  roll-position  stabilised, 
causes  rolling  until  the  demanded  manoeuvre  (e.g.  curvature  of  the  fli^t 
path  in  one  plane J  is  in  a  plane  inclined  at  45°  to  the  wing  planes. 

If  there  is  marked  rolling  of  the  missile  during  homing,  miss  distances 
may  be  appreciably  increased.  Tests  indicated  roll-position  stabilisation  to 
be  preferable  to  stabilisation  of  roll  rate  in  such  cases. 

Aerodynamic  non-linearities  and  cross-couplings  appear  to  have  little 
adverse  effect  on  stability,  provided  that  the  missile  has  a  suitable  auto¬ 
pilot.  Cruciform  missiles  so  provided  should  be  capable  of  manoeuvre  at 
incidences  up  to  30°  at  least. 

The  autopilot  of  V.R.725  appeared  to  be  satisfactory  from  this  point  of 
view,  despite  some  evidence  from  flight  trials  that,  at  least,  transitory 
instability  in  roll  may  occur  when  the  missile  incidence  exceeds  15°  in  a 
median  plane.  The  cause  of  such  behaviour  was  not  discovered,  but  may  have 
been  unsteady  air  flows  which  were  not  observed  or  simulated. 

In  the  configuration  studied,  the  non-linearities  in  pitch  and  yaw 
occurred  mainly  in  the  aerodynamic  moments,  and  affected  response,  particu¬ 
larly  at  incidences  exceeding  15°.  When  the  relative  v;ind  was  not  parallel 
to  a  wing  plane,  the  response  v/as  further  modified  by  aerodynamic  cross¬ 
coupling  between  the  pitch  and  yaw  planes,  and  by  the  non-linear  rolling 
moment  referred  to  above. 

A  constant  demand  for  a  slow  rate  of  roll  -  to  even  out  the  effects  of 
fortuitous  bias  -  appeared  to  be  detrimental  to  the  performance  of  a  missile 
free  from  bias,  miss  distances  being  increased  in  general. 

A  limited  investigation  of  the  effects  of  noise,  in  the  presence  of 
aerodynamic  characteristics  of  the  type  mentioned,  indicated  that  angular  noise 
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tended  to  increase  miss  distance  without  obscuring  the  effects  of  rolling 
which  have  been  mentioned  above. 

A  comparison  of  digital  and  analogue  solutions  showed  that  TRIDAC 
reproduced  the  response  and  stability  characteristics  of  the  mathematical 
model  with  generally  high,  and  always  adequate,  precision,  but  did  not 
generate  miss  distances  sufficiently  accurately  for  the  effects  of  aero¬ 
dynamic  non-linearities  on  homing  performance  to  be  established  precisely. 
The  effect  on  miss  distance  being  generally  small,  resort  to  digital  com¬ 
putation  was  necessary  to  obtain  the  differentieils  with  sufficient  accuracy 
for  definite  conclusions  to  be  drawn. 
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APPENDIX  1 
EQUATIONS  OF  MOTION 


Let  GjXYZ  (Fig, 2)  be  a  set  of  principal  axes  of  inertia  fixed  in  the 
missile,  the  origin  of  the  system  being  at  G,  the  centre  of  gravity.  Then, 
if,  in  these  axes, 

u,v,w  are  components  of  linear  velocity  of  the  missile 
p,q,r  are  components  of  angular  velocity  of  the  missile 
X,Y,Z  are  components  of  force  on  the  missile 
L,M,N  are  components  of  moment  on  the  missile 
A,B, C  are  the  moments  of  inertia, 

and  m  is  the  mass  of  the  missile,  assumed  constant,  at  any  rate  over  a  short 
interval  of  time,  the  equations  of  motion  are 


u  + 

qw  -  rv  =  X/m 

(1) 

V  + 

ru  -  pw  =  Y/m 

(2) 

w  + 

pv  -  qu  =  Z/m 

(3) 

Ap  - 

(B-C)qr  =  L 

(4) 

B4  - 

(C-A)rp  =  M 

(5) 

Cr  - 

(A-B)pq  =  N  . 

(6) 

If  V  is  the  total  velocity  of  the  missile,  so  that 

Tr2  2  2  2  /_v 

V  =  u  +  V  +  w  ,  {7) 


it  follows,  on  multiplying  equations  (1),  (2)  and  (3)  by  u,  v,  and  w, 
respectively,  and  adding,  that 


W  =  1  (Xu+Yv+Zw) 


(8) 


whence,  if  V  is  constant. 


Xu  +  Yv  +  Zw  =  0  •  (9) 

Either  equation  (8),  or  its  particular  case  (9),  may  replace  equation  (1)  and 
be  used  to  define  the  X  force,  including  the  motor  thrust,  aerodynamic  drag 
and  gravity  component,  necessary  to  maintain  any  given  variation  of  V  with 
time,  including,  as  a  special  case,  V  =  constant.  This  was  done  in  the 
simulation.  In  the  case  of  a  symmetrical  cruciform  missile,  it  may  be 
assumed  that  B  =  C.  If,  then. 


k  = 


the  equations  (A) »  (5)  and  (6)  reduce  to 


(10) 
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(11) 

q  -  krp  =  M/B  (12) 

r  +  kqp  =  N/C  (13) 

These,  with  equations  (2)  and  (3),  and  V  =  f(t),  where  the  function  f(t)  is 
defined,  were  the  equations  used  to  specify  the  motion.  The  forces  and 
moments  were  assumed  to  'be  wholly  aerodynamic,  gravity  components  being 
neglected.  For  the  most  part,  they  were  assumed  to  be  completely  defined  by 
the  static  aerodynamic  characteristics,  as  measured  in  the  wind  tunnel. 

The  inertia  terms  arising  from  the  motions  of  the  controls  and  com¬ 
ponents  of  the  control  system  were  neglected,  i.e.,  the  body  was  assumed  to 
be  rigid,  although  the  dependence  of  the  forces  and  moments  on  the  control 
positions  was  included. 

The  angular  positions  of  the  controls  were  related  to  the  angular 
velocities  of  the  missile  and  the  demanded  rates  through  the  transfer  func¬ 
tions  of  the  control  system. 

Axis  transformations 


p  =  L/A 


The  transfoimation  of  the  target  motion  from  earth  to  missile  axes  was 
accomplished  by  setting  up  a  three  gimbal  system  in  which  the  orders  of  rota¬ 
tion  were  pitch,  yaw,  roll,  and  the  corresponding  angles  of  rotation  were 
Yg^,  Zg  and  Xg,  The  relationships  between  the  co-ordinates  in  the  two  systems 

are  given  in  Ref.lZf  (Appendix  III,  Section  2,1,  Configuration  l),  and  those 
between  the  missile  and  gimbal  angular  velocities  are  also  quoted  (Section  3» 
Configuration  1 ) . 

The  transf onnation  of  the  components  of  relative  velocity,  and  of  the 
missile  spins,  from  missile  to  dish  axes  was  accomplished  by  setting  up  a  two 
gimbal  system  in  which  the  order  of  rotation  was  yaw  followed  by  pitch,  and 
the  corresponding  angles  of  rotation  were  0^^.  The  resulting  relationships 

between  the  components  in  the  two  systems  are  the  same  as  those  between  co¬ 
ordinates  quoted  in  Ref. 14  (Section  6,4,  Homing  Equations,  or  Appendix  'VI, 
equations  (VI , 2) ) . 

The  spins,  and  of  the  dish  about  its  gimbals  are  equivalent  to 
spins  -  ijfjj  sin  0^^,  0^^,  ijfjj  cos  0^^  about  its  principal  axes.  Hence,  the  total 
spins  of  the  dish  about  its  axes  are: 

Pg  =  (p  cos  iJfjj  +  q  sin  \lfjj)  cos  0^^  -  (r+i^Tj^)  sin  0^^ 

q^  =  -p  sin  tjj  +  q  cos  +  0^^ 

Tg  =  (p  cos  ijfjj  +  q  sin  sin  0jj  +  (r+Jj^)  cos  0^^  , 

These  spins  were  employed  in  the  integration  of  the  relative  velocities 
of  target  and  missile  in  dish  axes  into  relative  positions,  in  accordance  with 


the  standard  equations : 

^  = 

*H  *  -  '■.yR 

^R  = 

h  *  “•3=^  - 

II 

4  *  Ve  - 1.=^ 

where  the  subscript  R  denotes  "relative". 
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APPENDIX  2 

THE  (t>  SERVO 


As  shown  in  Fig. 2,  the  angle  is  defined  by  the  relation  tan  <f>  =  v/w, 
and  the  requirement  was  for  a  servo  to  generate  the  angle  (j)  continuously 
when  supplied  with  the  variables  v  and  w.  The  resolvers  on  the  shaft  would 
then  generate  the  sines  and  cosines  of  (j)  and  4*^.  The  problem  is  the  well- 
known  one  of  Cartesian  to  polar  co-ordinate  transformation,  and  the  method 
is  a  standai'd  technique  although  it  had  not  previously  been  used  on  TRIDAC. 

By  using  one  sine  and  one  cosine  channel  on  the  4>  shaft  (to  which  the 
h4>  shaft  was  ganged  mechanically),  a  servo  error  signal  e  is  generated 
from  the  relation 


e  =  V  cos  (^  -  w  sin  (j>. 


(14) 


In  the  steady  state,  when  e  =  0,  this  gives  the  required  condition 
tan  (j>  =  v/w.  The  error  is  brought  to  zero  by  arranging  that  the  shaft 
velocity  is  proportional  to  e  and  in  a  direction  to  reduce  it. 

To  do  this,  the  servo  used,  and  represented  by  K.&(p).  in  Fig.32(a), 
was  a  standard  TRIDAC  hydraulic  closed-loop  velocity  servo'^3.  The  overall 
loop  is  of  the  type  often  termed  "implied  position",  and  the  arrangement 
differs  from  the  normal  position  servo  in  being  highly  non-linear  in 
operation. 

Using  further  sine  and  cosine  channels  on  the  <f>  shaft,  the  servo 
provided  the  quantity  V0  from  the  relation: 

V0  =  V  sin  0  +  w  cos  <^.  (15) 


When  V  was  variable,  0  was  obtained  from  V0  by  division^^  using  a 
further  seiTVo;  but,  in  most  of  the  work,  V  was  constant  and  0  was  obtained 
merely  by  scaling  V0  suitably. 

From  (14),  it  can  be  seen  that,  when  v  and  w  are  both  small,  the  error 
signal  will  be  small  even  v,hen  (f>  is  substantially  in  error,  and  may,  in  fact, be 
insufficient  to  drive  the  servo.  It  is  also  apparent  that  (14)  is  satisfied 
in  equilibrium  by  a  value  of  (f>  differing  by  18O®  from  the  required  value. 

This  would,  from  (15),  produce  negative  values  of  5.  As  a  servo  problem, 
these  difficulties  can  be  illustrated  by  a  rearrangement  of  the  servo  loop. 


Define  a  servo  input  angle 

sin  d),  = 

cos  (f>.  = 


such  that 

v/^v  -fW 

1/  2  2  \  ^ 
w/^v  +w  } 


(16) 


where  the  square  root  has  the  positive  sign. 

The  servo  output  angle  being  it  is  required  to  make  follow  <l>^. 
Substituting  (16)  in  (14),  it  follows  that: 
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=  (v^  +  sin  -  <f)^) 

=  (v^  +  sin  (f)  (17) 

£ 

^  d>,  (b  • 

If,  now,  the  error  drives  a  closed-loop  velocity  servo  of  transfer 
function  K/p  (at  low  frequencies),  i.e., 

“^o  _  K 

e  ~  p  " 

we  have  ^ 

6  =  —  (v^  +  sin  6  •  (18) 

Hence,  the  re-arranged  loop  is  as  in  Pig, 32(b). 

Pig*32(c)  is  a  plot  of  the  error  function  of  equation  (17)  and  the 

normal  linear  position-servo  error-function,  shown  by  the  dotted  line  AOA. 

For  stability  equal  to  that  of  a  linear  position  servo  formed  by  feedback  at 

fixed  gain  around  the  velocity  servo,  the  slope  at  the  oapigin  of  the  sinusoid 

corresponding  to  the  maximum  possible  value  of  (v2^+  w2)'2  must  be  that  of  the 

linear  servo  characteristic.  Thus,  when  (v^  +  w2)'§'  is  less  than  its  maximum, 

and/or  (f)  is  not  small,  the  effective  slope  is  reduced  and  performance 
£ 

suffers.  The  point  B  is  obviously  one  of  unstable  equilibrium,  since  the 
slope  of  the  characteristic  is  negative  there.  Any  disturbance  from  B  will 
produce  a  velocity  towards  0,  either  clockwise  or  anticlockwise,  according  to 
the  direction  of  the  disturbance.  Failure  to  do  this  would  result  in  a  nega¬ 
tive  0  as  well  as  an  error  of  180°  in  (^  ,  and  hence,  <f)  . 

In  the  above  equations,  the  quantities  v  and  w  are  voltages  representing 
the  velocities  of  sideslip,  which  combine  to  replace  the  normal  single  feed¬ 
back  voltage  (derived  from  a  fixed  reference  voltage)  of  a  linear  position 
servo.  If  the  full  scale  values  of  v  and  w  are  represented  by  the  sam^ 

voltage  as  the  above  reference  voltage,  the  maximum  value  of  (v^  +  w^)^  is  V2 

times  that  voltage.  Thus,  under  conditions  of  maximum  input  and  small  servo 

error,  the  loop  gain  would  be  V’2  times  that  of  the  linear  servo,  and  gain 

adjustment  at  some  other  point  in  the  loop  is  required  if  the  same  degree  of 
stability  is  to  be  maintained. 

To  combat  the  falling  loop  gain  for  small  inputs,  a  form  of  automatic 
gain  control  was  fitted  to  the  servo,  whereby^ the  gain  of  one  amplifier  in 
the  loop  was  raised  as  the  value  of  (v2  +  wty^  fell._  IVhen  the  shaft  is  in 
the  correipt  position,  the  servo  forms  the  quantity  V0  which  equals 
(v^  +  w^)'2.  A  small  auxiliary  electric  servo,  therefore,  was  driven  by  this 
voltage  and  arranted  to  control  the  main  loop  gain,  as  indicated  in  Fig. 32(d) . 
Since  rapid  changes  of  v  and  w  were  not  encountered  in  the  problem,  the 
relatively  poor  dynamic  performance  pf  this  small  servo  was  immaterial.  Full 
compensation  for  values  of  (v2  +  w2)'2  right  down  to  zero  was  not,  of  course, 
possible,  nor  vms  it  necessary.  A  range  of  gain  of  20:1  was  found  to  be 
about  the  limit  in  practice.  Nevertheless,  <f)  was  placed  in  the  right  quad¬ 
rant  provided  that  the  value  of  (v2  +  w^)'2  equalled,  or  exceeded,  0*3  per 
cent  of  full  scale.  At  the  minimum  value,  the  angle  <j)  was  correct  within 


where 
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half  a  degree.  The  effect  of  this  error  on  the  computed  value  of  6^ is 
extremely  small  and  could  not  be  measured.  For  values  of  (v^  +  below 

0*3  per  cent  full  scale,  the  angle  0  could  not  be  held  at  or  near  its 
correct  value,  the  servo  tending  to  drift  at  random.  Under  these  conditions, 
however,  6  is  sensibly  zero  and  the  servo  had  no  function  to  perform. 

The  above  compensation  does  not  alter  the  theoretical  possibility  of 
unstable  point  operation.  Thus,  for  example,  in  single  plane  working  where 
V  may  change  sign,  if  the  angle  cf)  is  correctly  defined  for  one  particular 
sign,  the  duty  of  the  servo  is  to  switch  as_rapidly  as  possible  through  180° 
as  V  passes  through  zero,  thus  maintaining  6  positive.  In  theory,  there  is 
no  signal  to  do  this,  since  e  is  always  zero  but,  in  practice,  there  are 
small  drift  voltages  within  the  loops  which,  with  the  A.&.C.  system,  receive 
sufficient  gain  to  switch  the  servo  over, 

A  simple  and  yet  accurate  statement  of  the  servo's  dynamic  performance 
cannot  be  given,  realistic  inputs  being  very  difficult  to  specify  and 
reproduce,  A  test  was  made  in  which  v  and  w  were  sine  and  cosine  signals  of 
the  same  frequency  and  amplitude.  Under  these  conditions,  the  servo  shaft 
rotates  at  constant  speed  equal  to  the  circular  frequency  of  the  inputs,  and 
lags  behind  the  input  vector  by  a  fixed  angle  which  is  a  measure  of  performance. 
Unfortunately,  although  interesting  as  a  demonstration,  this  provides  no 
information  vjhich  is  not  already  known  from  the  constants  of  the  closed-loop 
velocity  servo  part  of  the  system. 

To  approach  more  nearly  an  operating  condition  found  in  the  roll 
stability  investigations,  the  v  and  w  channels  were  fed  with  equal  large 

D.C.  pedestal  voltages  to  which  small  sine-wave  signals  of  equal  amplitude, 

but  in  exact  anti-phase,  were  added.  Such  combined  inputs  produce  approximate 
simple  harmonic  motion  of  the  0  shaft  about  the  0  =  45°  position,  and  the 

frequency  response  method  was  employed  to  determine  the  lag.  For  inputs  con¬ 

sisting  of  1  c.p.s.  waveforms  of  r.m.s.  amplitude  0*5  per  cent  full  scale  on 
pedestal  voltages  of  50  per  cent  full  scale,  the  measurements  indicated  a 
rough  correspondence  to  a  simple  lag  of  about  2*D  m.s.  At  half  this  frequency 
and  at  an  r.m.s.  amplitude  of  1*2  per  cent  full  scale,  the  equivalent  lag 
was  roughly  halved. 

It  was  demonstrated  by  various  means,  such  as  the  incorporation  of 
phase  advance  external  to  the  servo,  increasing  loop  gain  by  a  factor  of  2 
and  decreasing  loop  gain  by  a  factor  of  10,  that  such  lags  as  were  present 
in  the  <p  servo  did  not  affect  the  overall  accuracy  of  the  simulation  signific¬ 
antly.  The  servo  v/as,  therefore,  considered  satisfactory. 
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SD.1ULAT0R  FAULTS  AND  ACCURACY 


The  fully  extended  model  required  some  400  amplifiers  to  simulate  the 
missile  homing  in  three  dimensions,  and  a  further  100  on  additional  require¬ 
ments  such  as  the  computation  of  miss  distance  and  total  lateral  accelera¬ 
tion,  and  for  noise  generation.  Most  of  the  available  capacity  and  special 
features  of  TRIDAC  were  required,  and  the  problem  of  keeping  the  computer 
"on  the  rails"  became  formidable.  Much  time,  therefore,  had  to  be  spent  in 
tracing  and  eradicating  faults  in  the  equipment. 

Two  main  types  of  fault  were  encountered.  The  more  difficult  to  cure 
was  the  intermittent  fault,  which  would  be  indicated  by  a  lack  of  repeat¬ 
ability  in  homing  runs  without  noise.  Such  faults  were  often  due  to  poorly 
mating  uniters  in  the  amplifier  brick  units.  The  other  main  type  of  fault 
was  that  due  to  the  total  or  partial  failure  of  a  component.  In  this  case, 
the  output  could  be  seen  to  be  incorrect  T«hen  a  standard  cheok  run  was  made 
and  the  results  compared  with  those  of  a  digital  solution.  If  the  fault  was 
not  indicated  by  the  automatic  monitoring  system,  as  frequently  happened, 
the  simulation  had  to  be  broken  down  and  tested  in  sections  until  the  faulty 
unit  was  isolated.  The  tracing  of  these  faults  was  made  more  difficult  by 
the  many  closed  loops  in  the  system.  Here,  the  digital  solutions  were  of 
great  assistance. 

When  the  simulator  was  producing  repeatable  results,  which  compared 
reasonably  well  with  the  digital  solutions,  the  main  concern  was  the  lack 
of  equality  of  the  miss  distance  values  in  homing  against  targets  performing 
symmetrical  manoeuvres.  For  example,  a  turn  of  the  target  to  right  or  left 
in  the  horizontal  plane,  following  a  head-on  approach,  would  result  in 
different  miss  distances,  the  values  depending  on  the  direction  of  turn. 

Prolonged  efforts  were  made  to  reduce  this  asymmetry,  but  the  final 
scatter  in  miss  distance  values  for  noise-free  targets  making  2g  tiirns  in 
either  horizontal  or  vertical,  planes  was  of  the  same  order  of  magnitude  as 
the  miss  distance  itself.  However,  within  the  limits  of  a  2g  target 
manoeuvre,  it  was  not  possible  to  obtain  miss  distances  greater  than  16  ft 
from  a  20,000  ft  range,  unless  noise  was  introduced. 

Fig, 30  compares  the  TRIDAC  results  with  those  obtained  digitally  in 
the  case  of  single  plane  homing.  The  initial  range  was  varied  from  20,000  ft 
downwards.  It  v;ill  be  seen  that,  in  spite  of  the  scatter  of  the  simulator 
results,  the  general  trend  of  the  digital  results  is  confirmed.  Pour  con¬ 
figurations  were  investigated,  namely,  target  turning  upwards,  or  downwards, 
in  the  vertical  plane,  or  to  port  or  starboard  horizontally.  Hence,  four 
values  of  miss  distance  were  obtained  for  each  initial  range.  At  values  of 
the  latter  exceeding  10,000  ft,  all  four  results  were  valid  but,  at  shorter 
ranges,  some  amplifiers  in  the  system  were  overloaded..  Re-scaling  to  over¬ 
come  this  was  not  considered  worthwhile  as  the  overloading  occurred  only 
when  the  amplifier  output  voltage  was  negative.  Consequently,  at  ranges  less 
than  10,000  ft,  there  are  only  two  values  of  miss  distance. 

The  digital  and  analogue  results  in  three  dimensions  are  compared  in 
Fig. 31  •  In  this  case,  the  initial  range  was  kept  constant,  but  the  initial 
roll  angle  was  varied.  At  each  value  of  the  latter,  four  values  of  miss 
distance  were  obtained,  namely,  two  by  simulating  the  manoeuvre  and  its 
mirror  image  in  the  yaw  plane,  and  two  by  repeating  the  process  in  the  pitch 
plane.  The  scatter  of  the  results  is  attributable  to  small  inaccuracies  in 
the  axis  transformations  viiioh  could  not  be  removed.  It  will  be  noted, 
however,  that,  over  much  of  the  range,  the  mean  of  the  four  values  at  each 
angle  is  close  to  the  digital  result. 
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The  miss  distances  shown  in  Figs, 30  and  31  were  measured  in  dish 
axes,  so  the  digital  results  shown  are  not  the  same  as  those  in  Figs, 23 
and  25,  which  represent  true  minimum  ranges. 

The  criterion  used  to  measure  the  asynmetries  of  the  simulation  was 
the  mEiximum  ohange  in  the  total  incidence  6  (expressed  as  a  percentage  of 
its  full  soale  value  of  0*6  radians)  due  to  changing  the  sign  of  the  input 
to  the  system.  In  the  homing  simulation  this  meant  comparing  a  symmetrical 
pair  of  target  manoeuvres.  In  the  ease  of  the  missile  alone,  it  meant 
changing  the  sign  of  the  demands.  When  biases  were  deliberately  introduoed, 
both  signs  of  bias  were  used  in  turn,  and  the  change_ln  5  was  divided  by 
twice  the  modulus  of  the  bias,  givir,g  the  ohange  in  9  per  unit  bias. 

The  asymmetries  are  shown  in  Table  3.  It  will  be  seen  that  the 
representation  of  the  missile  control  system  and  its  aerodynamios  was 
extremely  symmetrical.  When  the  homing  system  was  added,  and  a  stationary 
target  was  used,  so  that  the  earth  to  missile  transformation  was  not 
required,  the  asymmetry  in  0,  in  single  plane,  approached  0*5  per  cent, 
f.s,,  and  that  between  single  planes  at  90°  averaged  1*0  per  cent,  f.s. 

In  single  plane  homing  on  to  a  2g  turning  target,  the  asymmetry  was 
greater,  andhad  a  value  between  0*4  per  cent  and  1*7  per  cent,  f.s.,  depending 
on  the  plane  and  the  sense  of  the  turn.  Also,  the  difference  between  the 
two  single  planes  at  90°  Increased,  and  asymmetries  up  to  3  per  cent,  f.s. 
occurred.  With  homing  in  three  dimensions,  the  asymmetry  in  one  plane  was 
between  1*0  and  2*5  per  oent  f.s.,  and,  between  planes  at  90°  asymmetries 
of  up  to  8  per  cent,  f.s.  occurred.  The  repeatability  of  Individuail  results 
in  eO-l  runs  was  better  than  1*5  per  cent  f.s. 

It  was  suspected  that  much  of  the  asymmetry  was  due  to  the  resolution 
of  the  sine  and  oosine  potentiometers  on  the  hydraulic  servo-resolvers, 
which  is  limited  to  one  part  in  900.  To  investigate  this,  two  series  of 
tests  were  oonducted  after  all  the  sine  and  oosine  potentiometers  had  been 
carefully  set  up  and  trimmed,  and  the  gains  and  time  constants  of  the  system 
checked.  The  first  tests  were  concerned  with  the  intrinsic  symmetry  of  the 
resolvers  and  the  effects  of  the  limited  potentiometer  resolution,  while  the 
second  set  investigated  the  effects  of  small  biases  deliberately  introduced 
into  the  simulation. 

The  symmetry  of  the  servo-resolvers  was  ascertained  by  reversing  the 
sign  of  the  drive  to  eaoh  resolver  in  turn  and,  after  correcting  the  signs 
of  the  outputs  from  the  sine  potentiometers,  running  the  single  plane  homing 
with  the  particular  servo-resolver  in  the  reversed  condition.  As  shown  in 
Table  4,  in  three  out  of  the  four  cases,  the  asymmetry  was  affected  by  up  to 
1  per  oent  full  soale. 

The  terms  of  the  axis  transformations  which  were  thought  to  be 

j,,  in  the  earth  to  missile 

transformation,  and  u,^^  sin  iltp  and  cos  sin  0^^  in  the  missile  to 

dish  axis  transformation.  These  terms  were  suspeot  because  appreciable 
voltages  were  normally  applied  to  the  ends  of  the  sine  potentiometers,  and 
the  wire  to  wire  movement  of  the  wiper  at  small  shaft  angles  produced 
noticeable  steps  in  the  output  voltage.  Also,  setting  up  the  potentiometer 
on  the  resolver  unit  oould  only  be  done  to  within  ±  one  half  txirn,  owing  to 
the  position  of  the  centre  tap. 

In  order  to  assess  the  effect  of  these  wire  to  wire  steps,  the  sine 
terms  were  approximated  to  by  setting  sin  Yj,  =  Y^,,  etc,  and  the  multiplica¬ 
tions  were  accomplished  on  the  improved  electric  servo-multipliers'^"' 
referred  to  elsewhere  in  this  Note,  The  approximation  made  it  possible  to 


inaccurate  were  u^^  sin  Y^,  and  u,^^  cos  "5^  sin  Z 
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improve  the  accuracy  by  increasing  the  range  of  angular  displacement  of  the 
servo  shaft,  thus  effectually  reducing  the  size  of  the  vdre  to  wire  steps  by 
increasing  their  number.  (The  resolution  of  the  servo  is  0*1  per  cent.)  An 
improvement  in  symmetry  of  some  0*75  per  cent  (full  scale)  res\ilted  but, 
since  large  gimbal  angles  occurred  during  homing,  the  approximation  could 
not  be  retained  for  general  use. 

An  attempt  was  then  made  to  improve  matters  by  using  two  sine  potentio¬ 
meters  in  pai’allel.  This  effected  a  noticeable  improvement  in  symmetry  in  the 
case  of  sin  and  the  additional  potentiometer  was  installed 

permanently.  A  further  improvement  was  made  by  replacing  the  wire-wound 
feedback  potentiometers  on  two  of  the  servo-resolvers  by  carbon  film 
potentiometers  of  the  same  basic  linearity  but  of  much  hi^er  resolution. 

This  was  done  in  the  case  of  both  servo-resolvers  used  in  the  pitch  plane  in 
single  plane  simulation  and  an  improvement  in  symmetry  of  0*5  per  cent  f.s. 
was  observed.  Unfortunately,  no  more  of  these  potentiometers  were  available 
at  the  time. 

A  static  check  of  the  complete  earth  to  missile  axis  transformation 
was  made  by  setting  the  gimbal  angles  at  preselected  values  and  applying  a 
suitable  input  voltage  in  one  co-ordinate  at  a  time.  The  output  voltages  in 
all  cases  were  accurate  to  within  0*5  per  cent  of  the  full  scale  value 
(i.e.  the  input  voltage).  This  compares  with  a  static  accuracy  of  0*3  per 
cent  f.s.  for  individual  resolver  units  tested  by  Inspection  Depairtment, 
R.A.E. 


Biases  were  introduced  into  the  simulation  of  single  plane  homing  in 
both  single  planes,  and  also  when  simulating  homing  in  three  dimensions. 

The  biases  were  intended  to  represent  possible  errors  in  setting  up  the  sine 
and  cosine  potentiometers,  and  errors  in  the  alignment  of  the  feedback 
potentiometers  of  the  servo-resolvers.  Each  of  these  errors  could  amount  to 
one  half  turn  of  the  wire  on  the  potentiometer.  The  effects  of  these  biases 
are  presented  in  terms  of  the  maximum  change  in  the  value  of  6  per  turn  error 
in  the  setting  of  the  potentiometer,  expressed  as  a  percentage  of  the  full 
scale  value  of  5,  when  homing  on  a  target  making  a  2g  turn  from  a  head-on 
approach.  A  further  type  of  error  simulated  by  the  biases  was  drift  in  the 
voltage  representing  the  angular  velocity  drive  to  the  seirvo-resolvers.  The 
effects  of  this  are  presented  as  the  maximum  change  in  0  per_millivolt  of 
drift,  expressed  as  a  percentage  of  the  f\ill  scale  value  of  6. 


The  significant  results  are  listed  in  Table  5*  These  show  that,  for 
two  dimensional  homing,  the  alignment  of  the  dish  servo  shaft,  and  of  its 
sine  potentiometer  operating  on  had  an  important  effect  on  the  results  - 

accuracy  being  important  since  u_„,  was  large  -  while  misalignment  of  the 

axis  transformation  servo  shaft  and  its  sine  potentiometer  operating  on 

had  a  lesser  effect.  Drift  in  the  drive  to  the  dish  servo  had  no  effect, 
while  a  small  change  was  caused  by  drift  in  the  servo  of  the  earth  to  missile 
transformation . 


In  three  dimensional  homing,  the  most  important  effect  was  due  to 

TEE 


misalignment  of  the  sine  potentiometers  forming  Urnrm  sin  and 


u_ 


TEE 


cos  Yj,  sin  Zj,  in  the  earth  to  missile  axis  transformation  of  the  target 


velocity.  This  affected  the  results  to  the  extent  of  causing  over  2  per  cent 
maximum  change  in  0,  expressed  as  a  percentage  of  full  scale  per  turn  of 
misalignment . 


The  next  most  important  effect  was  that  due  to  biasing  of  the  sine 
potentiometers  involved  in  generating  the  components  of  dish  spin  rates 
depending  on  the  missile  rodl  rate,  but  it  was  an  order  of  magnitude  smaller 
than  the  previous  one.  All  the  other  biases  which  were  applied  to  the 
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the  servo- resolver  drives  and  shaft  initial  positions,  the  sine  potentio¬ 
meters  of  the  missile  to  dish  axis  transformation,  the  earth  to  missile 
axis  transformation  of  v_rr,  and  w__,  and  the  generation  of  components  of 

p  -  the  dish  axial  spin  -  from  the  missile  pitch  and  yaw  rates,  had  a 
s 

minor  effect. 

These  tests  indicated  that  the  lack  of  symmetry  in  the  simulation  was 
mostly  due  to  inaccuracies  in  certain  terms  of  the  axis  transformations 
containing  sines  of  the  gimbal  angles.  A  secondary  cause  of  error  was 
probably  small  drifts  in  parts  of  the  model  not  in  a  feedback  loop.  The 
drive  to  the  earth  to  missile  axis  transformation  is  a  case  in  point. 

It  is  probable  that  the  accuracy  of  the  axis  transformations  could  be 
improved  by  the  use  of  carbon  film  pencil  potentiometers  in  critical 
positions.  A  prototype  is  under  development.  A  further  improvement  is  to 
be  expected  from  fitting  carbon  film  feedback  potentiometers  to  all  the 
hydraulic  seivo-re solvers  and  multipliers. 

More  recently,  it  has  come  light  that  drifts  have  been  occurring  in 
some  integrators,  which  have  not  been  indicated  by  the  monitoring  system. 
These  may  have  explained  some  day  to  day  variations  in  the  results  which 
caused  much  trouble  and  loss  of  time.  The  short  term  stability  of  the 
simulator  was  generally  very  good. 
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Basic  data  on  D.4  and  W.2  (V.R.725)  as  simulated 


Units 

D.4 

W.2 

(V.R.725) 

Missile  Dimensions  Length 

Feet 

18-8 

20*8 

Body  Max.  Diameter 

Wing  Span 

Feet 

Feet 

1-71 

5-33 

1-75 

5-33 

Distance  from  datum(wing  tip  leading 
edge) to  tail  hinge 

Feet 

5* 

5*96 

Reference  Area  S 

Feet^ 

1-767 

1-767 

Reference  Length  o 

Feet 

1-5 

1-5 

Mass 

At  end  of  boost 

Slugs 

60-7 

73-5 

At  end  of  sustainer 

Slugs 

39-8 

53-1 

Usod  in  simulation 

Slugs 

59-8 

59-9 

Moment  of  Inertia  A 

Used  in  simulation 

Slug  ft^ 

36*6 

31  -2 

Moments  of  Inertia  B  and  C 

Used  in  simulation 

Slug  ft^ 

887-1 

1177 

Distance  of  C.G-.  ahead  of  datum 

^wing  tip  leading  edge) 

Used  in  simulation 

Inches 

6-3 

8*  5 

Shift  during  burning  of  sustainer 
motor 

Inches 

1*9 

5-0 

Maximum  incidence  of  wind  tunnel  tests 

(Nominal)  M  =  1*6 

22° 

35° 

M  =  2-0 

25° 

40° 
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Control  surface  deflection  coefficients 


^  =  -v^ 

h  (e  <f>  <T^)  =  +D^  (6  .#.  +  90  cr^) 

\  (6  ^  <T^)  =  -D^  (0  .#.  +  90  cr^) 

\  (0  4>  o-^)  =  +D^  (0  ^  +  90  o-^) 

(0  (f>  (X^)  =  +D^  (0  (f)  +  180  a^) 

\  (0  cp  0-3)  =  -D^  (0  ^  +  180  cr^) 

(0  ^  0-3)  =  -D^  (0  ^  +  180  o-^) 

(0  cl>  0-^)  =  +D^  (0  ^  +  270  o-^) 

^  o'i,)  =  +^ni  ^  ^'l^ 

^n  ^0  ^  =  ”^n  ^'l^ 


(0  (j>  -cr^)  =  -D^  (0  18O  -  (f)  a^) 
\  (6  <#>  -o-^)  =  (0  180  -  <!>  o-^) 

\  (0  <#>  -^1)  =  +Dn  (0  180  -  ^  o-^) 

(0  "0'2^  ~  90  -  i#>  cr^) 

(0  <!>  ■<T^)  =  -D^  (0  90  -  o-^) 

(0  .#.  -<T2)  =  -D^  (0  90  -  cr^) 

(6  ^  ~  ^ 

^  ~  ■'■^  O'-]) 

Dn  (0  <#>  -O3)  =  ^'l^ 

(0  <!>  -o-^)  =  -D^  (0  270  -  <l>  a^) 

\  ^  ^'l^ 

Dn  (0  -o-^)  =  +D^  (0  270  -  .#,  a^) 


5  =  Uo-^  +  o-^)  =  Uo-2  + 

o-^  =  Tl  +  ? 

T)  =  -  o-^) 

0-2  =  S  +  ? 

S  =  ^(og  -  cr^) 

O-J  =  -T1  +  € 

0-^  =  -s  +  ? 

D^jDffljDn  represent  the  changes  in  the  respective  moment  coefficients 

resulting  from  the  deflection  from  zero  of  one  control  panel.  Couplings 
between  panels  are  neglected. 
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TABLE  3 

Ohserved  asymmetries  in  6,  maximum  change  in  6;  full  scale 

(39S  F.S.  t  1°  change  in  6) 


Target  Manoeuvre 

0 

Dimen- 
sional 
•  com¬ 
plexity 
of 

problem 

Asymmeti*ies 

In-Plane 

Between 

Planes 

at 

90° 

Missile  alone 

3 

0-0*1 

O'*  1-0*3 

Offset  stationary  target  (range  = 

20  K  ft) 

2  . 

0*3-0* 5 

0*5-1 *5 

2g  target  1  " 

"  } 

2 

0*4-1 *7 

2*0-3*0 

2g  target  (  " 

"  ) 

3 

1 *0-2*5 

3*0-8*0 

Repeatability  with  2g  target 

3 

0-1*5 

0-1*5 

TABLE  4 

Effect  of  modifications  on  asymmetry,  maximum  change  in  6;  ^  full  scale 


Modification 

Dimensional 

Complexity 

of 

Problem 

Maximum 

Change 

Notes 

Reversal  of  dish  servo  in 

(a)  Yaw, 

2 

0*5 

(b)  Pitch 

2 

1*0 

Reversal  of  gimbal  servo  in 

(a)  Yaw, 

2 

0 

(b)  Pitch 

2 

0*5 

Replacement  of  sine  potentio¬ 
meters  by  high  performance 
electric  servo  multipliers  in 
fa)  Yaw, 

2 

0*75 

Appr oximati on, 
only  valid  for 
small  shaft 
angles 

(b)  Pitch 

2 

0*75 

Paralleling  sine  potentiometers 
on  dish  servo  resolver-'in 
(a)  Yaw, 

2 

2*0 

Permanently 

installed 

(b)  Pitch 

2 

0 

Paralleling  sine  potentiometers 
on  gimbal  servo  resolver  in 

fa)  Yaw, 

2 

0 

(b)  Pitch 

2 

0*5 

Carbon  film  feedback  potentio¬ 
meters  fitted  to  the  servo 
resolvers  in  Pitch  plane 

2 

0*5 

Outputs 

considerably 

smoother 
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TABLE  5 

Asymmetries  due  to  applied  Liases 

Maximum  change  in  6;  %  full  scale 


Quantity  Biased 

Dimensional 
Complexity 
of  Problem 

In-Plane  Asymmetry 
Resulting 

Uj^j,  sin  (gimbal  angle) 

2 

0-2  per  turn 

Gimbal  position  (initial) 

2 

0*2  per  turn 

Gimbal  drive 

2 

0*2  per  mV 

'^TMM  angle) 

2 

0*6  per  turn 

Dish  initial  position 

2 

0*6  per  turn 

Dish  drive 

2 

0  per  mV 

Utee  (gimbal  angle) 

3 

1 •4-2*2  per  turn 

p  sin  (dish  angle) 

3 

0* 1  per  turn 

\ 
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FIG.  I. 


FIG.  I.  APPROXIMATE  HALF  SECTIONS  OF  D4  &  VR725. 
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FIG.2.  DIAGRAMMATIC  REPRESENTATION  OF  MISSIl  F 
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FIG.  I.  APPROXIMATE  HALF  SECTIONS  OF  D4  &  VR725. 
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FIG.  I. 
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FIG.  3. 


FIG.  3.  POLYNOMIAL  FITS  TO  D-4 


PITCHING  MOMENT  DATA. 

MACH  No.  1-58 
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FIG.4.(a&hO 
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FIG.  4.  (a. «  AUTOPILOTS  USED  IN  THE  SIMULATION. 
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FIG.  5. 
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FIG.  5.  BLOCK  DIAGRAM  OF  D4  SIMULATION. 
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FIG.  6. 
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TAIL- FIXED  AERODYNAMICS  OF  VR.  725  MODEL. 


FIG.6.  F-FUNCTIONS  CONTRIBUTING 
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FIG.7.  F- FUNCTION  CONTRIBUTING  TO  cf,  . 
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FIG.8.  F-FUNCTION  DEFINING 
INDUCED  ROLLING  MOMENT. 
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FIG. 9. 


F -FUNCTIONS  CONTRIBUTING 
TO  NORMAL  FORCES 
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VARIATION  OF  D- FUNCTIONS 
WITH  ROLL  ANGLE  <p. 
6=20°  M- 1-6,  a-IO°20°30°. 

(see  table  2) 
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FIG.  II.  BLOCK  DIAGRAM  OF  VR  725  SIMULATION  AT  CONSTANT  MACH  No., 
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FIG.  12  RESPONSE  OF  MODEL  IN  SINGLE  PLANE  TO  RATE  OF  PITCH  DEMANI 

WITH  LINEAR  AND  NON-LINEAR  AERODYNAMICS. 
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FIG.  13  YAW  PLANE  RESPONSE  OF  MODEL 
TO  EQUAL  YAW  &  PITCH  RATE  DEMANDS. 
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FIG.  15.  RESPONSE  OF  MODEL  IN  3  DIMENSIONS  TO  3  DEMANDS  FOR  ANGULAR  VELOCITY 

IN  ROLL  STABILISED  AXES. 

NON-LINEAR  BODY-WING  AND  LINEAR  CONTROL  SURFACE  AERODYNAMICS. 
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FIG.  16.  RESPONSE  OF  MODEL  IN  3  DIMENSIONS  TO  3  DEMANDS  FOR  ANGULAR  VELOCITY 

IN  ROLL  STABILISED  AXES. 

NON-LINEAR  BODY -WING  AND  CONTROL  SURFACE  AERODYNAMICS. 
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FIG.  17  COMPARISON  OF  SIMULATION  WITH 
TELEMETRY  RECORDS  FOR  ROUND  035. 
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FIG. 18. 


FIG.  18.  DETAILS  OF  CONTROL  ROUNDS  INCLUDING  TELEMETERED  ROLL  RESPONSE. 
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FIG.  19  &  20. 


SIMULATED  ALTITUDE  40  K  FT. 
SIMULATED  MACH  No.  1-6 
SIMULATED  INCIDENCE  E7-5® 


FIG.  19.  DEPENDENCE  OF  LAG  AND  ROLL 
OSCILLATION  FREQUENCY  ON  L^  . 
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FIG.  20.  EFFECT  OF  LAGS  ON  CRITICAL  INCIDENCE. 
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FIG.  21 
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FIG  21.  BLOCK  DIAGRAM  OF  HOMING  MISSILE  SIMULATION. 
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FIG. 22.  BLOCK  DIAGRAM  OF  HOMING  HEAD  IN  ONE  PLANE. 
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FIG.  2  3.  a  2  4. 


FIG  23.  MISS  DISTANCE  v.  INITIAL 
RANGE  IN  SINGLE  PLANE  HOMING. 


FIG.24,  PEAK  LATERAL  ACCELERATION  v. 
INITIAL  RANGE  IN  SINGLE  PLANE  HOMING. 
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FIG.  25.(a&b) 


INITIAL  ROLL  ANOLE  X^q 


(a)  INITIAL  RANGE  14,000  FT. 


(b)  INITIAL  RANGE  20,000  FT. 


FIG.  25.(a*b)  MISS  DISTANCE  v.  INITIAL  ROLL 
ANGLE  IN  3-DIMENSIONXli:  HOMING  FROM  TWO 

INITIAL  RANGES. 
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FIG.26.  EFFECT  ON  MISS  DISTANCE 

OF  A  STEADY  ROLL  RATE 


FIG.  27.  EFFECT  ON  MISS  DISTANCE 

OF  ANGULAR  NOISE 
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FIG.29.  RESPONSE  IN  3  DIMENSIONS: 
COMPARISON  OF  DIGITAL  AND  ANALOGUE  RESULTS 
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FIG.30.  COMPARISON  OF  ANALOGUE 

AND  DIGITAL  RESULTS  IN  SINGLE 
PLANE  HOMING. 
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FIG  31.  COMPARISOM'-OF  ANALOGUE  AND 
DIGITAL  RESULTS  IN  THREE  DIMENSIONS. 
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(a)  GENERATION  OF  0  AND  6 


(b)  RE-ARRANGED  SERVO  LOOP. 


(c)  ERROR  CHARACTERISTIC 


FIG.32Ca-d) SERVO  ARRANGEMENT  FOR 

GENERATION  OF  0  AND  0 
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